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ABSTRACT 

An analysis of the flow across a sealing dam of the type that appears in gas turbine 
seals was developed for steady, laminar, subsonic, isothermal compressible flow. The 
analysis is valid for both parallel sealing dam surfaces and surfaces separated by a 
small tilt angle. The computer program determines mass flow rate, pressure and 
velocity distributions, Mach number, force, center of pressure, axial film stiffness, 
rotational flow and pressure flow Reynolds numbers, power loss, and approximate 
temperature rise resulting from viscous shearing. The output is in both English units 
and the In L er national System of Units. Some of the results can be automatically plotted. 
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COMPUTER PROGRAM FOR VISCOUS, ISOTHERMAL COMPRESSIBLE 
FLOW ACROSS A SEALING DAM WITH SMALL TILT ANGLE 
by John Zuk and Patricia J. Smith 
Lewis Research Center 

SUMMARY 

An analysis of the flow across a sealing dam of the type that appears in gas turbine 
seals was developed for steady, laminar, subsonic, isothermal compressible flow. 

Both parallel sealing dam surfaces and surfaces separated by a small tilt angle (e. g. , 
due to thermal distortion) are treated. The analysis is valid for plane pressure flow. 
When the sealing dam mean radius is much greater than the sealing dam radial width 
which, in turn, is much greater than the mean film thickness, the analysis is valid for 
hydrostatic radial flow. The analysis is also valid for relative rotation of the sealing 
dam surfaces if the circumferential shear flow has little effect on the radial pressure 
flow. 

A computer program to carry out this analysis is also presented. Input variables 
include the dimensions of the seal, pressure boundary conditions, and molecular weight 
and physical properties of the gas. The output includes mass flow rate, pressure and 
velocity distributions, Mach number, force, axial film stiffness, center of pressure, 
rotational and pressure flow Reynolds numbers, Knudsen number, torque, power loss, 
and approximate temperature rise resulting from the viscous shearing for specified film 
thicknesses. The output units are in both the English and International Systems. Some 
of the results can be automatically plotted. 


INTRODUCTION 

Some powerplants, such as advanced jet engines, exceed the operating limits for 
which face contact seals were designed (refs. 1 and 2). As a result, noncontact face 
seals have become necessary. A noncontact face seal which is pressure (force) balanced 
is shown in figure 1. In this seal, the pressure drop occurs across a narrowly spaced 
sealing dam, and the force due to this pressure drop is balanced by a predetermined hy- 






Figure 1. - Pressure-balanced face seal with no axial film stiffness. 


drostatic closing force and a spring force. This configuration, however, has an inherent 
problem. For parallel sealing dam surfaces, there is a force caused by the pressure drop 
across the sealing dam, and this force is independent of film thickness; hence, there is 
no way of maintaining a preselected film thickness which will allow tolerable leakage and 
still have noncontact operation. Since the force is independent of film thickness, the 
design also lacks axial film stiffness for sufficient dynamic tracking of the stationary 
nosepiece with the rotating seal seat. The seal nosepiece must follow the seal seat sur - 
face under many conditions without surface contact or excessive increase in film thick- 
ness, which would yield high leakage. Some of these conditions are axial runout, mis- 
alinement, thermal distortion, coning, and dishing. 

A promising method of maintaining a preselected film thickness and achieving axial 
film stiffness is to add a gas bearing, such as a shrouded Rayleigh step pad bearing, to 
the noncontact pressure -balanced seal (refs. 1 and 2). This is illustrated in figure 2. 

Both the sealing dam force due to the pressure drop across the sealing dam and the gas 
bearing force are balanced by the hydrostatic and spring closing forces. The gas bearing 
has a desirable characteristic whereby the force increases with decreasing film thick- 
ness. If the seal is perturbed in such a way as to decrease the gap, the additional force 
generated by the gas bearing will open the gap to the original equilibrium position. In a 
similar manner, if the gap becomes larger, the gas bearing force decreases, and the 
closing force will cause the seal gap to return to the equilibrium position. 

Since a proper balance of the opening and closing forces must be found in order to 
determine a gap with a tolerable mass leakage, physical quantities of interest such as 
pressure distribution and mass leakage must be calculated for several gaps. The pres- 
sure distribution and mass leakage have been calculated for the parallel film hydrostatic 
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Figure 2. - Pressure-balanced face seal with a gas bearing added for axial film stiffness. 

case. Mathematical solutions for the hydrostatic, isothermal, compressible, viscous 
flow sealing dam exist in the literature (e. g. , see Gross (ref. 3)). Carothers (ref. 4) 
has conducted compressible flow experiments in thin films of air which are in radial flow 
between parallel plates. The pressure distribution was found for both subsonic and super- 
sonic axial entrance flows. Grinnell (ref. 5) has theoretically and experimentally investi- 
gated compressible flow in a thin passage and has shown excellent agreement between 
theory and experiment. 

An analysis for small tilts of the sealing dam surfaces does not exist in the litera- 
ture. The tilts are caused by mechanical and thermal distortions of the seal, as stated 
in reference 1. To achieve a good design it is desirable to study the effect of the variation 
of a large number of parameters; thus, an automatic calculation and printout of physical 
variables would facilitate design. 

This report presents a compressible flow sealing dam analysis that considers both 
parallel sealing dam surfaces and surfaces with small tilt angles. A computer program 
was developed to perform the calculations for the sealing dam design. Input variables 
include seal dimensions, pressure boundary conditions, and molecular weight and physi- 
cal properties of the gas. Physical quantities such as mass flow rate, pressure and 
velocity distributions, Mach number, force, power, torque, and center of pressure are 
found for specified film thicknesses. 
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BASIC MODEL AND EQUATIONS 


The sealing dam model consists of two parallel, concentric, circular rings in rela- 
tive rotation at a constant speed separated by a very narrow gap. A pressure differen- 
tial exists between the rings' inner and outer radii (see fig. 3). 

The model formulation is based on the following physical conditions: 

(1) The fluid is homogeneous, compressible, viscous, and Newtonian. 

(2) The flow is steady and laminar (continuum flow regime), and the body forces are 
negligible. 

(3) The bulk modulus is ignored (X = -2/3 p). This is Stokes idealization (ref. 3). 
This condition will be valid unless the gas is under high pressure, very dense (e. g. , 
shock wave structure), or rarefied. 

(4) The fluid behaves as a perfect gas. 

(5) Since AR is much greater than h, the entrance region effects are neglected; 
hence, the convective inertia forces are neglected. This means that the seal is treated as 
operating entirely in the viscous region. 

(6) The fluid film is isothermal. This means that all heat generated in the film is 
conducted away through the walls. This is a standard assumption of lubrication theory. 



Figure 3. - Model of the sealing dam with a small tilt angle (not to scale). 
Rotating upper ring removed for clarity. 
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The validity of this assumption breaks down for cases of large thermal gradients in the 
sealing dam and when the frictional heating is high (e. g. , small gap or high speed). 
However, thermal analysis of a seal (unpublished data by T. E. Russell of Lewis) shows 
that the sealing dam can be closely approximated by a constant temperature. In any case, 
the Mach number must be less than 1 /-yfy. This is the limit of the validity of isothermal 
duct flow analyses as stated in most gas dynamics textbooks (e. g. , Shapiro (ref. 6)). 

(7) The entrance Mach number is close to zero. This means that there cannot be a 
large axial -flow source on one surface impinging on the radial surface, as is present in 
a hydrostatic bearing. 

(8) The fluid velocity in the reservoir is considered to be negligible (stagnant) and thus 
its effects are neglected in this analysis. 

(9) The radial pressure flow is uncoupled from the rotational shear flow. Reference 7 
shows this assumption is valid under the following conditions: 

(a) The ratio of the reference rotational velocity to the reference radial flow 
velocity must be less than 1/ ^Re^h/AR) (i. e. , Rfi/U ref < l/\^Re h (h/AR)). 

(b) The mean radius R must be much greater than h, and Re^h/AR) must be 
much less than 1. 

Using the above two conditions, reference 7 further shows 

(c) The treatment of the radial flow as uncoupled from the rotational flow is a 
very good approximation for most applications where the radial pressure differential 
is large and the speeds are moderate. 

(d) The circumferential and axial pressure variations are negligible. 

(3) If AR/R^ is much less than 1, a two-dimensional channel or narrow slot is 
valid. 

The rotational flow, however, is important for power loss calculations due to viscous 
shear and transition to turbulent flow. 

The governing flow equations for a compressible fluid with constant viscosity in vec- 
tor notation are (ref. 8) 

Conservation of mass: 


5£ + p V • V = 0 
Dt 


Conservation of momentum (Navier -Stokes equations): 


p — = - VP - /i Curl (Curl V) + (X + 2p) V (V • V) + F 
Dt 

Isothermal equation of state: 


P = P(p) 
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All symbols are defined in appendix A. 

A rectilinear Cartesian coordinate system is used to describe the radial sealing dam 
flow (see fig. 3). 

Applying the conditions assumed in the model reduces the above system qf equations 
to the following set which can be solved: 

Conservation of mass: 

3 (£!!) = 0 (1) 

0X 

This form of the continuity equation is not used but is replaced by the integrated 
form, which is shown to be the conservation of mass flow in the radial direction. 

Conservation of momentum: 


3 ^ 

— -=0 z-direction (circumferential) 


— = q ^ x -direction (radial) 

* ay 2 

Equation of state: 

P Pi 

p = p <ft T = A 

Pi 

Solving the circumferential direction momentum equation (eq. (2)) gives 


(2) 

(3) 


(4) 


v = Cj(r)y + C 2 (r) 


Applying the boundary conditions 

v= 0~ C 2 = 0 at y = 0 

v=rf2-C 1 = — at y=h 

1 h 


yields 
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v = 

h 


(5) 


Since the azimuthal and axial pressure variations are neglected, P is a function of x only 
(P = P(x)), and the radial momentum equation can be integrated twice with respect to y 


1 dP 2 ^ ^ t 

u = y + C.y + Co 

2/i dx 1 * 


Applying the boundary conditions 

u = 0 — Cg = 0 at y = 0 

u.0-cS-®£ at y = h(x) 
1 2 dx 

yields 

u = — — (y 2 -hy) 

2 ii dx 

Now, the mass flow at any x per unit width is 


M 

L 



h 3 (x)p dP 
12 p dx 


Substituting the perfect gas law (eq. (4)) into equation (7) yields 

M _ h3p l dP 2 
L 24^ dx 


( 6 ) 


(7) 


( 8 ) 


The mass flow does not change with x; hence, 

— = 0 (9) 

dx 

(This equation replaces the continuity equation (1).) The boundary conditions for pressure 
are 

P = P , or P 2 = P 2 
P = P 2 , or P 2 = P 2 

An alternate formulation of this problem would have been to start with the compres- 
sible Reynolds’ lubrication equation (ref. 3) for this model. Reynolds’ equation is of the 
form of equation (9) with boundary conditions (10). 


at x = 0 


at x = R 2 - R x 


( 10 ) 
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Parallel Film Limiting Case 


When the seal surfaces are parallel, h is constant, and equation (9) is readily solved 
by integrating twice and applying the boundary conditions (eq. (10)). The result is 



The mass leakage is found by substituting equation (12) into equation (8). This yields 


M = 


Lh 3 p 1 P 1 
24p(R2 — k 




(13) 


Equations (11) and (13) are reducible to the form found in reference 3. The total force 
per unit width is found from 


F 

L 


/ 


R 2” R 1 


(P - P - )dx 
v min' 


which results in 


F 

L 



(14) 
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Where P m j n is the smaller pressure of the two pressure boundary conditions. The 
center of pressure is 


L / R2 ' R 1 


(P - P min )x dx 


X c = 


( 15 ) 


which results in 



Since equation (14) does not depend on film thickness, there is no axial film stiffness for 
parallel sealing surfaces. This deficiency is overcome by adding a gas bearing to the 
sealing faces (fig. 2). 


Small -Tilt-Angle Case, a 4 0 

^or nonparallel sealing dam surfaces, the film thickness is no longer a constant; 
thus, equation (8) must be solved for the variable -film -thickness case. For a sealing 
dam with a small tilt (see fig. 4), the mean clearance is 
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For small tilts, 


sin a = 


h(x) - h 
x - C 


m 


= a 


where 


R 0 - R 


C = 


2 iV l 


Hence, h(x) = h m + a(x - C) or h(x) = B + ax where B = h m - a C. If 


D = - 


24jxP i 


Pi 1 


equation (8) becomes 


DM 


dP 


* (ax + B) 3 


Now the conservation of mass must be satisfied, that is, dM/dx = 0 or 


d 2 P 2 -3aDM 


dx 2 (ax + B) 4 


(17) 


Integrating equation (17) twice yields 


p 3 = -_ DM + c '; x + c ' 2 ' 
2a(ax + B) 2 


(18) 


When the following boundary conditions are applied, 
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P = P x - C 2 = P 2 + at X = 0 


P = P 2 -Ci = 


p2 p2 

*2 *1 


2aB 


DMh 


m 


l 2" R l B 2 (h m + aC) 




at x = 2C 


( 19 ) 


equation (18) becomes 


p 2 _ DM 

2a 


B 2 (ax + B) 2 


+ P 


1 - 


P 2> 

II 


Pj/\ r 2 ' R iy 


DMh x 
m 


B 2 ( h m + aC) 2 


( 20 ) 


or 


p 2 = p 2 


i - i - 


2 

,2 


\ R 2 " j 


DM 


B 


(ax + B) 2 - B 2 
2a(ax + B) 2 


h x 
m 


<v + “ c ) 2 


(21) 


For small tilts of the sealing dam surfaces, the mass leakage is found from equa- 
tion (13), where h is now that is, 


Lh 3 , p 1 P 1 / P 3> 

M = char 1 1 [i 2 

24/i(R2 “ Rj) l p2 


(22) 


where 


char 


h l h 2 


m 


Using equation (22) for the mass leakage in equation (21), the pressure distribution for 
small tilt angles becomes 


\ 
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( 23 ) 



24 mM (ax + B) 2 - B 2 
PjPjLB 2 2a(ax + B) 2 


The total force per unit width is found from 


h x 
m 


1/2 


(h m ^c) ; 


F 

L 



(P - P . )dx 
v min' 


(24) 


The center of pressure is found from 


X 


c 



-R, 


(P - 


P . )xdx 
min' 


F 

L 


(25) 


In the computer program, the integrations in equations (24) and (25) are performed 
numerically., 

The axial film stiffness is defined as 


STIFF = 


dF 


dh 


m 


(26) 


It can be seen from equation (23) that a sealing dam with nonparallel surfaces has an 
axial film stiffness. It can also be seen that this film stiffness can be positive or nega- 
tive. In the computer program, the axial film stiffness is found numerically (see appen- 
dix B). 


Additional Parameters Calculated by Computer Program 

The average radial velocity at any radial point x in the sealing dam gap is found 
from 


u (x) = — — — 
av p(x)hL 


MP 1 

Lp 1 P(x)h 


(27) 
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The local Mach number at any x is then 


M = 


u av< x ) 


u av< x > 

VyfflT 


(28) 


where a is the speed of sound. 

The pressure flow Reynolds number is found by using the hydraulic radius 2h as the 
characteristic length 


Re h = 


P 2 u av 2h 

ju(RT 


The Knudsen number can be found from (ref. 9) 


(29) 


1 48 M 

_ Molecular me an free path ^ j max 

Mean film thickness Re^ 

2 


(30) 


Under conditions of very small film thicknesses, the Knudsen number may be greater 
than 0. 01, and this continuum analysis would no longer be valid. (A slip flow regime 
analysis must be used. ) 

The total power is found by considering only the viscous shear caused by rotation. 


- r 2 o 2 

Power = RQ • (Shear force) = — - — ^ 

h 


/ 


dA 


_ jiR 2 S2 2 A 


(31) 


where RS2 is the mean rotational velocity. 

A very rough estimate of the film temperature rise due to the viscous shearing can be 
estimated by equating the heat generated by viscous shearing with heat transfer by convec- 
tion. Thus, 


I 


rr _ ^ _ M A 2 A 

film, av 

hC p M 


(32) 
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This calculated film temperature rise will be higher than actual film temperatures. 
The predominant mode of heat transfer, conduction by the walls, is neglected. 


Computer Program Formulation 

The previously derived equations, placed in the form used in the computer program, 
are shown in table I in the English system of units. A description of the program and flow 
charts is presented in appendix B. The program listing is given in appendix C. A sam- 
ple problem with its input and output for parallel sealing dam surfaces and relative sur- 
face tilts of ±1 milliradians is given in appendix D. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 9, 1969, 

120-27-04-90-22. 
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TABLE I. - FORM OF PERTINENT EQUATIONS FOUND IN COMPUTER PROGRAM (IN SAME SEQUENCE AS EQUATIONS APPEAR IN PROGRAM) 
A = n (r 2 - R 2 ), in. 2 


<R = — , 
M 


(lbm)( R) 
P 


- where (R = 1545, 4 


ft -lb 


(Ib-mole)(°R) 


Pi" 




<R(T + 460) 


32. 174 


(lbm)(ft) 

(lb)(sec 2 ) 


( lb)(sec z ) 

ft 4 


l/y(R(T + 460) 

32 174 

f 

(lb)(sec 2 ) 


ft 

sec 


R-» + Rq 

If L = 0 , L = 2 tt -± 1 , in. 


“ p 0 = °> Po' Pi 


(l b)(sec^) 
' ft 4 


If N * 0, V = ( R - 1 — ( min \ iL 

[ 12 in. 1 \ 2 J \ 60 sec/’ sec 


If N = 0 and V = 0, 
AT = 0, °R 

Btu 


C = 0, 


Re(R) 


(Ibm)(°R) 




/60 sec\ /144 inA 

\ min / \ ft 2 J 


Number of steps 


, in. 


m = jit 


M = 32. 174 


( Ibm)(ft) 

(Ib)(sec 2 ) 



h char p l P l L j 1 — 
l, 


24p(R 2 - Rj) 


min 


h min “ smaller of h j and h 2 > in * 


MP. 


32 174 0k m )(ft) 

/60 sec\l 

(lb)(sec 2 ) 

\ min / 


ft 

P(x) ’ sec 


M = - 


Re(P) = 


[32. 174 ( lb m )( ft ) "j m r(t + 4 
L (Ib)(sec 2 )J 


Q = 13. 083 M, std cu ft/min 
Power 


\ h m / SSOjNIb 12_irJ 


[42. 42 ( lbXmin)l ( p ower ) 

L b P J 


H t0(al = 42.42(Power), 

min 


/*Rr>-R 

F = L * HP ~ P^dx, lb 

T /"Rq-R.. 

x c -7 J '( p - p 0 > x <*• «"• 

F "0 

x = x 1 + n Ax, in. 
h = h m + “( x ' RA 7^ ' in - 

— =1/ t - L - — ^ /■ X in. 

Pi r \ pyb'i 


If a = 0, -£-=1 


If a * 0, — = < 



DM 

(ox + B) 2 -_B^ 

h x 
m 

} 

y p 2 y \ r 2 ■ R i/ 

B 2 p2 

_ 2n(ox + B) 2 

(h m + c,C)\ 

1 


dF lb^ 
dli > in. 


|P t -P 2 |(R 2 -Ri)L 
X 


Torque = J 33 0Q0 -C I) g^_ er ) (l* jj^ , ft- 
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APPENDIX A 


A 

a 

B 

C 

C 

C 

P 

C v 

D 

F 

F 

F 

g 

h 

K 

Kn 

L 

M 

M 

m 

n 

P 

P min 

Q 

R 

R 

AR 


SYMBOLS 

2 2 

sealing dam surface area, in. ; m 
speed of sound, ft/sec; m/sec 
h m -aC 

constant of integration 

(R 2 - R^/2 

specific heat at constant pressure, Btu/(lb)(°R); J/(kg)(K) 
specific heat at constant volume, Btu/(lb)(°R); J/(kg)(K) 
-24pP 1 /p 1 L 

sealing dam force, lb; N 
dimensionless force, F^Pg - P^)(R 2 - R-^) 
body force vector 

gravitational constant, 32.174 (lbm)(ft)/(lb)(sec^); 9.81 m/sec^ 
film thickness, in. ; m 
DM 

Knudsen number 
sealing dam width, in. ; m 
Mach number 

mass flow, lb/min; kg/sec 

molecular weight of gas, lbm/lb-mole; kg/kg -mole 
an integer 

2 

static pressure, psi; N/m 

2 

smaller pressure of two pressure boundary conditions, psi; N/m 
net volume flow rate, std cu ft/min; std cu m/sec 
radius, in. ; m 

mean radius, (R^ + R^/Z, in. ; m 
sealing dam length, Rg - R^, in. ; m 
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(R gas constant, universal gas constant/molecular weight, ft-lb/(lbm)(°R); 

J/(kg)(K) 

(R universal gas constant, 1545. 4 ft-lb/(lb-mole)(°R); 8. 3143 J/(kg-mole)(K) 

Re Reynolds number 

r radial direction coordinate 

STIFF axial film stiffness = -dF/dh m , lb/in.; N/m 
T temperature, °F; K 

U pressure flow reference velocity, ft/sec; m/sec 

u velocity in r -direction or x -direction, ft/sec; m/sec 

V moving sealing dam surface speed (circumferential direction), ft/sec; m/sec 

V fluid velocity vector 

v velocity in circumferential direction or z-direction, ft/sec; m/sec 

X c center of pressure in radial or x-direction, in. ; m 

X dimensionless center of pressure, X /(R 0 - R.) 

x coordinate in pressure gradient direction (x-direction) 

y coordinate across film thickness 

z shear flow coordinate in Cartesian system 

a relative inclination angle of sealing dam surfaces, rad 

y specific heat ratio, C /C 

r ^ 

9 circumferential coordinate or azimuthal direction 

X second viscosity coefficient or coefficient of bulk viscosity 

2 2 

jU absolute or dynamic viscosity, (lb-sec)/ft ; N -sec/m 

p density, (lb)(sec 2 )/ft^; kg/m 3 

Cl angular rotation velocity, rad/sec 

V Del operator, — * + — ^ 

3x 3y 3z 

Subscripts: 
av average 

char characteristic 

h based on film thickness 

m mean 
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max 


maximum 


min minimum 

r based on radius 

0 reference 

1 inner radius or inlet 

2 outer radius or outlet 
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APPENDIX B 


COMPUTER PROGRAM 

The program called SEAL performs an analysis of the flow across a gas film sealing 
dam with mean film thickness h m and tilt angle a. SEAL and its subprograms are 
written in FORTRAN IV. (The computer at the Lewis Research Center is an IBM 
709411/7044 or 7040 Direct Couple computer under IBSYS version 13 using ALTIO. ) 

Included in this appendix are lists of the input variables and program variables, 
along with their descriptions and English units; detailed descriptions of the subprograms; 
and a flow chart of the main program (fig. 5). 

All input, calculations, and output are in English units. Printout of data in Inter- 
national Units is optional. 


Main Program 

The main program, SEAL, performs the primary flow analysis. Subroutines are 
used for secondary operations such as numerical integrations, numerical differentiation, 
and plotting data. 

Input to SEAL is by punched cards in the following order: 

(1) Title card - alphanumeric identification of the data (format 12A6) 

(2) NJ card - number of mean film thicknesses to be analyzed in one running of the 

program (format 13) 

(3) h m cards - mean film thicknesses, six per card (format F12. 6) 

(4) Data cards - seal dimensions, pressure boundary conditions, physical properties 

of the gas, and logical variables (read by NAMELIST/INPT/) 

Data are read by NAMELIST to minimize the number of cards required to run a second 
case with the same title and h m cards. Input variables are initially set to zero. Con- 
sequently, variables which are not changed during the reading of the data cards will be 
calculated by the program. (See the list of input variables for their significance and for 
any restrictions on them. ) 

Output data are printed in English units in groups, in the following order: 

(1) Program identification - compressible sealing dam with small tilt angle 

(2) Data identification - as it appears on the title card 
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Do statements 200 to 400 for J «= 1, N J 



(a) Initial steps. 


(b) Main calculation. 


Figure 5. - Flow chart of main program (SEAL). 
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(c) Write routine. 


°; 0 > — — — ~@> 


No 



(d) Plot routine. 


Figure 5. - Concluded. 
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(3) Input data - as it appears on INPT cards 

(4) Calculated constants - rotational flow Reynolds number, reference density, 

density at inner radius of seal, speed of sound, gas constant, length of seal, 
and rotational velocity 

(5) Parameters that vary with h m - mass flow rate, volume flow rate at standard 

conditions, maximum Mach number, pressure flow Reynolds number, rotational 
flow Reynolds number, sealing dam force, axial film stiffness, center of pres- 
sure, dimensionless sealing dam force, dimensionless center of pressure, and 
Knudsen number (Printout of dimensionless sealing dam force and dimensionless 
center of pressure may be suppressed by setting RSKIP to TRUE. ) 

(6) Parameters associated with power dissipation - power, shear heat, apparent 

temperature rise, and torque (Printout of all parameters associated with power 
dissipation may be suppressed by setting TSKIP to TRUE. ) 

(7) Parameters that vary across the sealing dam (one group for each h ) - distance 

across sealing dam, pressure, pressure ratio, average velocity, and Mach 
number (Printout of all data in group 7 may be suppressed by setting ASKIP to 
TRUE.) 

Printout of data in English units is followed by plots of several parameters (see 
appendix D). Plots appear in standard form with minimum x and minimum y in the 
lower left corner of the plot. Legends at the bottom of the plots give conversion factors 
for International Units. The plots appear in the following order: 

(1) Center of pressure as a function of h m 

(2) Dimensionless sealing dam force as a function of h m (Plots 1 and 2 are sup- 

pressed for a = 0 since center of pressure and sealing dam force do not change 
with h . For a * 0, plot 2 may be suppressed by setting RSKIP to TRUE. ) 

(3) Power as a function of h m for nonzero speed (Plot 3 may be suppressed by setting 

TSKIP to TRUE. ) 

(4) Pressure ratio as a function of distance across sealing dam for subsonic flow 

cases (If a = 0, only one plot is made because the pressure distributions are 
identical for all h . All plots in group 4 may be suppressed by setting ASKIP 
to TRUE. ) 

Following the plots, data are printed in International Units. This printout may be sup- 
pressed by setting NOUI to TRUE. 

SEAL is divided roughly into seven sections. The first section reads data and calcu- 
lates program constants (cards 41 to 59). 
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The second section (cards 85 to 93) tests certain input variables. If they are zero, 
new values are calculated. If they are nonzero, the original values are used by the pro- 
gram. Since SPEED and CAPV both represent rotational velocity, they must be consist- 
ent. If SPEED is read as zero, CAPV must be examined. In the case that CAPV is not 
zero, SPEED is calculated from CAPV. If they are both zero, the system is considered 
to be static. 

Section three (cards 99 to 146) is the first part of a loop which is done for each h . 
This section calculates starting values for x and h, mass flow rate, volume flow rate 
at standard conditions, pressure flow Reynolds number, rotational flow Reynolds number, 
maximum Mach number, and Knudsen number. If the maximum Mach number indicates 
that the flow analysis is no longer valid (M > 1/Vr), IHTAG(J) is set equal to 1 as a trig- 
ger, and no further calculations are made. If the flow analysis remains valid, the pro- 
gram calculates power, shear heat, temperature rise due to power dissipation, torque, 
sealing dam force, and center of pressure. For a + 0, the integrations in the force and 
center of pressure equations are done numerically by Simpson’s rule. For a = 0, the 
integrations are done analytically, and the resulting formulas are used in the program. 
The axial film stiffness is found numerically by Lagrange differentiation of force with 
respect to h . 

Section four (cards 151 to 163) is the rest of the loop started in section three. Sec- 
tion four calculates film thickness, pressure, pressure ratio, average velocity, and 
Mach number at several points across the sealing dam. Section five (cards 168 to 211) 
writes data in English units. Section six (cards 215 to 257) plots the various parameters. 
And section seven (cards 261 to 319) writes data in International Units. 

Numerical constants that appear in the program are for units conversion. 


Subprograms 

SEAL uses six subroutines whose listings are given in appendix C. They are SEMPS1, 
PX, PXX, PRESS, STFNSS, and ARRNG. The last five are described in table II. SEAL 
also uses two subroutines that are not standard in IBSYS. These are SORTXY and 
PLOTXY. 

SORTXY sorts two numerical arrays. A statement such as CALL SORTXY(X, Y, N) 
results in the array X being rearranged such that X(l) < X(2) . . . ^ X(N) with the 
Y array rearranged to preserve the X, Y pairs. N is the number of elements in the 
X and Y arrays. 
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TABLE n. - DESCRIPTION OF SUBROUTINES 


Name 

Call vector 

COMMON block 

Progran 

Description 


variables 

variables 

variables 


PRESS (pressure function) 

X 

PI 


distance from inner radius of seal 
pressure at inner radius of seal 



P2 


smaller of PI and P2 



alpha 


tilt angle 



c 


1 - P^/pS 



RDIF 


distance across seal 



HM 


mean film thickness 



AK 


constant 




YY 

distance from inner radius of seal 




Y 

pressure 




AA 

film thickness at outer radius of seal 




Q 

P 2 for a * 0 

PX (integrand in integral) 

Y 



distance from inner radius of seal 

r R 2 -r i 


PI 


pressure at inner radius of seal 

A < P - P mm )dx 


PREF 


smaller of PI and P2 



A, B, C, D, E 

YY 

dummy variables to fill common block 
distance from inner radius of seal 

PXX (integrand in integral) 

r R 2~ R i 

/ (P - P . )x dx 

x min 

Y 

PI 


distance from inner radius of seal 
pressure at inner radius of seal 


PREF 


smaller of PI and P2 



A, B, C, D, E 

YY 

dummy variables to fill common block 
distance from inner radius of seal 

STFNSS (numerical differentia- 

XX 



input array of independent variable 

tion for film stiffness) 

YY 



input array of dependent variable 


ITAG 



input array of numerical flags 





(ITAG(J) * 0 implies flow became 
supersonic and case should be elimi- 
nated from calculation. ) 


DDY 



returned array of film stiffness 





DDY = -d(YY) / d(XX) 


MAX 


X 

number of elements in XX, YY, and ITAG 
array of valid independent variables 




Y 

array of valid dependent variables 




MM 

number of valid points 




N 

number of points used in the numerical 





differentiation, N - 5 




K 

index of point at which differentiation 





is made 




1ST 

index of first point used in differen- 





tiation 




IN 

index of last point used in differen- 





tiation 




II. JJ 

indices of points used in differentiation 





(II can be used as either the row 
index or the column index since 
the matrix A is square. ) 




A(5, 5) 

matrix whose elements A( I, J) = 





X(II) -X(JJ) 




PI 

product (x k - 3C.)P.(x.) 




SI 

sum f A. 




S2 

c 

3 

>T 1 
x- 




P2 

n 

product | | (x^ - Xj) 

j- i 




KY 

integer value of derivative 




DY(50) 

floating point array of derivatives 

ARRNG (arranges arrays to 

X 



nput array of independent variable 

be plotted) 

Y 


: 

nput array of dependent variable 


XP 



sorted array of new independent variable 


YP 



sorted array of new dependent variable 


N 



lumber of elements in input arrays 


I 


T |i 

lumber of elements in sorted arrays 
emporary storage for sorting 



PLOTXY plots two numerical arrays. A statement such as CALL PLOTXY (X, Y, 
KODE, P) produces an on-line plot of X against Y with the point X(l), Y(l) in the 
upper left corner. The plot appears with X increasing down the page and Y increasing 
across the page. KODE indicates which plotting options are used. For example, 

KODE = 6 gives a plot with most of the grid lines suppressed, * as the plotting character, 
and the X and Y scales computed by the plotting routine. The array P contains in- 
formation needed by the plotting routine, such as the number of points to be plotted, the 
X and Y scales if the programmer computes them, and the frequency of grid lines in 
the X and Y directions. 

Special format statements are used to print plot titles and plot legends. A pair of 
statements such as 


WRITE (6, 1) 

1 FORMAT (2HPT, 10HPLOT TITLE) 

will print the title PLOT TITLE above the plot, A pair of statements such as 

WRITE (6, 2) 

2 FORMAT (2HPL, 11HPLOT LEGEND) 

will print the legend PLOT LEGEND immediately below the plot. 

Listing of SORTXY and PLOTXY can be found in reference 10 or can be obtained 
from the Instrument and Computing Division of the Lewis Research Center. 

SIMPS1 is a function subprogram used to perform a numerical integration by 
Simpson’s rule. A statement such as F = SIMPS1 (XO, XF, G, K) gives F as the defi- 
nite integral 


F = 



G(X)dX 


The integrand is evaluated at interior points by the external function G named in the cal- 
ling vector. The interval of integration, XO to XF, is not divided uniformly. More 
subdivisions are made in regions where the integrand is changing rapidly. If two succes- 

_5 

sive evaluations of the integral on a particular subinterval differ by more than 3x10 
times the value of the integrand, the subinterval is divided into two subintervals and the 
integration is repeated. In the integration requires more than 200 subintervals, the 
integer K is raised by 1 to indicate that the returned value of the integral is incorrect. 

PRESS is a function subprogram to evaluate the pressure at any distance X from 
the inner radius of the sealing dam. The distance appears in the calling vector. The 
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pressure differential equation is solved analytically, and the resulting formulas are used 
in the program. For a = 0, the formula used is equation (11). 


P= Pjl/l - 


1 - 


p 2 R 2 - 


(ID 


For a 0, the formula used is equation (23) 

r_ _ 


P = P x < 


1_,1__2W * 


p 2 1 \ r 2 " R u 


24juM 

Pl LB 2 


(ax + B) 2 - B 2 


h m x 


~\2 


2o'(ofX + Bp (h m + aC) 


1/2 


(23) 


PX is an external function subprogram to evaluate the integrand in the integral 


F 

L 


/ 


R 2 - r i 


(P - P . )dx 
v min' 


(24) 


Similarly, PXX is an external function subprogram to evaluate the integrand in the 
integral 


/ 

•S r\ 


,R 2" R 1 


(P - P . )x dx 
v mm' 


X = 
c 


(25) 


F 

L 


STFNSS is a subroutine subprogram to calculate the axial film stiffness (STIFF). 
It performs the numerical differentiation 


STIFF = 

by Lagrange’s method (ref. 11). The general formula for Lagrange differentiation is 

n 

L’(x) = 

i=0 
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where 


and 



P.(x) = (x 




For the special case of the derivative at one of the tabulated points x^, 


L ' (x k> = JJ 

j*k 




1 



(34) 


where = (x k - x i )P.(x.), i * k. 

Subroutine STFNSS follows a computing scheme described in reference 11. A matrix 
is defined whose elements are 


a ij 


= x. 


- x. 
1 


The product of the off-diagonal elements of each row is multiplied by x, - x., except for 
th. ^ i 

the k row. This defines D^. The sum 



is formed and multiplied by the product of the negative of the off-diagonal elements in the 
j column. This gives the first term in the formula. The second term is formed by 
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■Hi 

summing the reciprocals of the off-diagonal elements in the k tn column and multiplying 
the sum by y^. 

The subroutine first eliminates cases for which the flow Mach number exceeds 
1/Vy. Then it arranges the data in order of ascending x. It chooses the five data points 
x i’ yi (i = k " 2, k - 1, k, k + 1, k + 2) for use in equation (34). If there are less than 
five data points in the set, all are used in the differentiation. 

ARRNG is a subroutine subprogram used to arrange two arrays X and Y for plot- 
ting. The subroutine first eliminates cases for which the Mach number exceeds 1/ Vy. 

It sorts the remaining data in order of ascending X. It then inverts both arrays and 
interchanges them. The data are now in the form of ordered pairs (Y(N),X(N)), (Y(N-l), 
X(N-l)), . . . , (Y(1),X(1)) where Y(N) > Y(N-l) > . . . > Y(l). Arranging data in this 
form permits the Y array to be plotted as the independent variable in descending order. 
Consequently, the plots appear in standard form with minimum X and minimum Y in 
the lower left corner. 


Input Variables 

Input variables to the program and their units are listed. Arrays are given with 
their dimensions. 


FORTRAN 

Unit 

Description 

symbol 

TITLE (12) 


alphanumeric identification of data 

NJ 


number of mean film thicknesses (NJ < 50) 

HMEAN(J) 

in. 

mean film thicknesses 

J=l, NJ 

ALPHA 

rad 

tilt angle 


L 

in. 

width of mean circumference of sealing dam 

SPEED 

rpm 

rotational speed'j 

Only one of these needs to 
be set. If both are set, 

CAPV 

ft/sec 

sealing dam 
surface speedy 

^ ' 
CAPV is calculated from 

SPEED. 

MOLWT 

lb/lb -mole 

molecular weight of gas 

PI 

psi 

pressure at inner radius of seal 
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FORTRAN 

symbol 

Unit 

Description 

P2 

psi 

pressure at outer radius of seal 

T 

°F 

isothermal reference temperature 

R1 

in. 

inner radius of seal 

R2 

in. 

outer radius of seal 

RHORO 

(lb)(sec 2 )/ft 4 

reference density at inner radius of seal (If 
RHORO is read as zero, program calculates 
RHORO. ) 

RHORF 

(lb)(sec 2 )/ft 4 

density at mean radius used in calculating ro- 
tational Reynolds number (If RHORF is read 
as zero, program calculates RHORF. ) 

MU 

(lb)(sec)/ft 2 

absolute viscosity of gas 

CP 

Btu/(lb)(°R) 

specific heat of gas 

GAMMA 


ratio of specific heats 

NGRID 


number of steps across seal (maximum, 20) 

ASKIP 


logical variable (If ASKIP = TRUE, program 
skips calculation and printout of x, pressure, 
average velocity, Mach number, and pres- 
sure ratio. It also skips plotting of x against 
pressure ratio. ) 

RSKIP 


logical variable (If RSKIP = TRUE, program 
skips calculation and printout of dimensionless 
center of pressure and dimensionless force. 

It also skips plotting of FBAR. ) 

TSKIP 


logical variable (If TSKIP = TRUE, program 
skips calculation, printout, and plotting of 
variables associated with power. ) 

NOUI 


logical variable (If NOUI = TRUE, program 


makes no conversion to International Units. ) 
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Program Variables 


The variables used in the program are listed in the approximate order of their ap- 
pearance. Arrays are given with their dimensions. Variables marked with a * are 
printed as output data. 

FORTRAN Unit Description 

symbol 


PI 

RUN IV 
ZERO 

ft -lb 

(lb-mole)(°R) 

PP1 

psi 

PREF 

psi 

AAA 

rad 

MCUT 


NN 


JMOD 


*AREA 

in. 2 

RDIF 

in. 

C 


CC 


PDIF 

psi 

*R 

ft -lb 


(lbm)( 0 R) 

♦RHOl 

(lb)(sec 2 )/ft 4 

A 

ft/ sec 


jt = 3. 1415927 
universal gas constant 

input variables are equated to ZERO which is 
initially set equal to 0. 0 

value of PI in common block 

smaller of PI and P2 in common block 

value of ALPHA in common block 

point at which mathematical model becomes 
invalid 

number of grid points (1 < NN ^21) 

number of pressure distributions that will fit 
evenly on one page 

face surface area of seal 

distance across seal, equal to R2 - R1 

outside- to inside -pressure ratio, equal to 
P2/P1 

constant used in pressure calculation 
(CC = 1 - (P2 2 /P1 2 )) 

total pressure drop across seal 
gas constant 

calculated density at inner radius of seal 
speed of sound 


FORTRAN 

symbol 

Unit 

Description 

IHTAG(J) 


numerical flag: 

IHTAG(J) = 0 implies flow analysis is valid 
IHTAG(J) * 0 implies flow analysis is not 
valid 

HM 

in. 

current value of HMEAN (J) in common block 

DELX 

in. 

distance between successive grid points 

*X(J,I) 

in. 

distance from inner radius of seal 

X(J, 1) 

in Q 

first point in x distribution 

X2 

in. 

last point in x distribution 

H(J, I) 

in. 

film thickness at X(J, I) 

H(J, 1) 

in. 

film thickness at X(J, I) 

H2 

in. 

film thickness at X2 

HCHAR 

in. 

characteristic film thickness 

*MDOT(J) 

lb/min 

mass flow across seal 

UAV 

ft/sec 

average velocity at outer radius of seal 

*MACHMX(J) 


maximum Mach number for given HMEAN(J) 

*REP(J) 


pressure Reynolds number 

RHOREF 

(lb)(sec 2 )/ft 4 

calculated density at mean radius of seal 

*RER(J) 


rotational Reynolds number 

*Q(J) 

std cu ft/min 

volume flow rate at standard conditions 

*KN(J) 


Knudsen number 

*POWER(J) 

hp 

power dissipated by viscous shearing 

*DELTJ(J) 

°R 

apparent temperature rise due to power dissipa 
tion 

*H TOTAL (J) 

Btu/min 

shear heat of system 

* TORQUE (J) 

ft-lb 

torque 

AK1 

(psi) 2 

constant needed in pressure calculations 
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FORTRAN 

Unit 

Description 

symbol 

K, KK 


numerical flags which indicate whether or not 

the numerical integrations of (P - P . ) and 

' min' 

(P - P . )x are accurate 
min 

*F(J) 

lb 

sealing dam force 

*XC(J) 

in. 

center of pressure 

*P(J, I) 

psi 

pressure at X(I, J) 

*PRAT(J, I) 


ratio of P to P . at X(I, J) 

*UAVRG(J, I) 

ft/sec 

average velocity at X(I, J) 

*MACH(J, I) 


Mach number at X(I, J) 

*STIFF(J) 

Ib/in. 

axial film stiffness, equal to -dF/dh m 

*FBAR(J) 


dimensionless force 

*XCBAR(J) 


dimensionless center of pressure 

JJ 


counter for cases for which Mach number 
< 1/ yfy (If JJ = 1 modulo JMOD, the printer 
will skip to a new page. ) 

XPT, YPT, 


utility arrays used in sorting and plotting 

XPLOT, YPLOT 

PP 


contains information needed by plotting subrou- 
tine(See ref. 10 for details. ) 

KODE 


plotting code (See ref. 10 for details. ) 

NP 


number of points in a plot (1 < NP < 50) 

*UI 


output variables in International Units 

I 


distance across sealing dam index 

J 


mean film thickness index 
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APPENDIX C 


PROGRAM LISTING 

$ I B F T C SEAL 1 

2 

COMPRESSIBLE FLOW SEALING OAM ANALYSIS WITH SMALL TILT ANGLE 3 

4 

LOGICAL ASKIP,RSKIP,TSKLP,.NUUI 5 

REAL MOOT , MOL WT,. MU, MACH, M ACHMX , MCUT , KN , L 6 

DIMENSION XPLOT < 50 ) „YP LOT ( 50 ) , XPT l 5 0 ) ,YPT(50) , PP{ 61) , ZER0()17) , 7 

1 Ul(30),TITLE(12). 3 

DIMENSION F ( 50) ,.XC ( 50) f MDOT ( 50 ) ,.XC B AR ( 50 ) , FBAR ( 50 ) , Q150 ) ,! 9 

1 POWER! 50) , HI ( 50) , HTOT AL ( 50) ,,DELT J (50 1 ) , TORQUE (50) 4 HME AN (50 ) * 10 

2 M ACHMX (50) , REP < 50 ) , S T I FF ( 50 ). , I HTAG ,( 5 0 > , K N ( 50 > , RER < 50 ) 11 

DIMENSION H <50*21 , X( 50,21) t-P ( 50,2 1 ) ,UAVRG ( 50 *2 1 ) , PRAT < 50 i2 1 1 * 12 

1 MAC H ( 5 0 , 2 1 ): 13 

COMMON/ I NTGR'L/PPl , PREF , A A A , CC , RDIF , HM,,AK 1 14 

EXTERNAL PX,PKX 15 

NAMELIST / I NPT/ ALPHA, .L,.SP&ED,CAPV, MOL WT , Pi, P2,T,RHR2*RH0RG, 16 

1 RHORF,MU,CP,.GAMMA,NGRID,ASKIP,RSKIP,TSKIP,NOUI 17 

DATA PI , RUN IV/ 3.1 4 1592 7, 1545. 4/ 18 

EQUIVALENCE ( Z ERO ( 1). , MOL W T ) , ( ZERO < 6 ) ,. P 1 ) , ( Z ERO i 1 1 > *L ) , 19 

1 ( Z. E RO ( 2 ) , ALPHA) , (ZER0t7), P2), < Z ERO t 1 2 ) 4 T ) , 20 

2 ( ZERO! 3), SPEED) , (ZERO* 8).. Rl), (ZERO 113 )‘#MU), 21 

3 (ZERQ(4),RHGR0), (ZERO (9)., R2) , l ZE RO 1 14 ) , C AP V ) , 22 

4 ( ZERO < 5 )., RHGRF ) ♦ ( ZERO ( 1 0 ) ,.C P ) , ( ZERO ( 1 5 ) , GAMMA ) 23 

DO 90 I =. 1 , 1 5 24 

90 Z-ERO (I). - 0. 25 

26 

READ LNPUT DAT A , C ALCUL A TE PROGRAM CONSTANTS, AND WRITE INPUT 27 

CONDITIONS 28 

29 

DATA CARDS 30 

TITLE - DATA IDENTIFICATION - 1 CARD (FORMAT 12A6) 31 

32 

NJ - NUMBER OF FILM THICKNESSES (FORMAT 13) 33 

34 

HMEAN - MEAN FILM THICKNESSES - 6 PER CARD (FORMAT 6FI2.0) 35 

36 

$ I.NPT - SEAL DIMENSIONS, OPERATING CONDITIONS# PHYSICAL 37 

PROPERTIES OF GAS, LOGICAL VARIABLES 38 

(READ BY NAMELIST/INPT/ ) 39 

40 

READ (5,-3) TITLE 41 

READ ( 5 ,. 1 ) NJ 42 

READ (5, .2) (HMEAN(J) ,J = 1,NJ) 43 

100 WRITE (6,10V 44 

READ ( 5 , I NP T) 45 

PP1 =i PI 46 

PREF = AMI.N1 (P1^P2) 47 

AAA ^ ALPHA 48 

MCUT = 1 • /SORT (GAMMA ) 49 

NN = NGRID+l 50 

J MOD = 59/<4+;NN) 51 

AREA = PI*(R2**2-R1**2) 52 

RDIF=iR2-Ri 53 

C =P 2 /P 1 54 

CC=1.-C*C 55 

PDIF^ABSt^l-Pa) 56 

R = RLNIV/MOLWT 57 

R HO 1 = Pl/R/(T+460.)*4. 4756636 58 

A— SQRT (GAMMA*R*(T+460.)*32.174) 59 

WRITE (6,54) TITLE 60 
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WRITE (6^11) ALPHA, P2,P1*T ,.MU,M0LWT f GAMMA,R2fRl»L*RH0R0,RH0RFrNN, 61 

l SP6ED„CAF>V*CP,ASKIPV R5KltP,TSKIP 62 

IE (MOUU GO TO 110 63 

64 

CONVERT I IMP UT DATA TO INTERNATIONAL UNITS 65 

66 

UI ( 1 K =5 ALPHA 67 

U I ( 2 ) P2*6..8947572E3 68 

UI( 3) = P1*6.8947572E3 69 

UK4). = ( T + 460 • 1/1 ,8 70 

U I ( 5 1 - MU * 4 7.88 0*2 5 8 71 

U I ( 6 ) =. MOL WT 72 

UI( 71 = GAMMA 73 

Ul ( 8 1 = R2*?.54E-2 74 

UI (9) =, Rl*2.54E-2 75 

UK 10) =. L*2 . 54E-2 76 

UIlll) = RH0RO*517.2026 77 

UI( 12) = RH0RF*517.2026 78 

UI ( 131 = SPEED/60 • 79 

UI ( 14)=. CAPV*. 3048 80 

UK 15) = CP*.41865783E4 81 

82 

TEST LNPUT PARAMETERS AND DETERMINE OPERATING CONDITONS 83 

84 

110 IF (L.EQ-O.) L-PI*(R1+R2> 85 

IF(RHORQ.NE.O. ) RH0l=RH0R0 86 

I.F( SPEED. EQ.O. > GO TO 120 87 

C AP V =i P l,*$PEED 3 MRl + R2) / 720. 88 

GO TO 200 89 

120 LFICAPV.NE.O.) GO TO 121 90 

TSKIP =. .TRUE. 91 

CP=0. 92 

121 SPEED = CAPV*720./PI/(Rl+R2) 93 

C 94 

C 98 

200 DO 4C0 vK 1 » N J 99 

IHTAGU) = 0 100 

HM = HMEAN(J) 101 

DEL X = R'D I.F /FLOAT ( NGR I D ) 102 

IF < ALPHA. GE.O. ) GO TO 210 103 

XU, 1) ~ RDIF 104 

X2 =0. 105 

DEL X = -DEL X 106 

GO TO 220 107 

210 X(J, 1) = 0. 108 

X2 = RDIE 109 

220 H ( J t 1 ) =. HMBAN( J)+ALPHA*(X( J, 1 )~R'DIF/2. ). 110 

H2 = HMEAN( J i )+ALPHA*(X2-RDIF/2.) Ill 

HCHAR = H(J , 1)*H2/BMEANU L 112 

MOOT U) = H;:HAR**3*RH01*P1*L*CC/MU/RDIF*6.7029 16 7 113 

UAV= ABSIMDOKJ ).) *P 1 /L /RH01 / AM I N1 IH ( J *1 ) ,.H2)/PREF/1 3 *405833 114 

MACHMX(J)= UAV/A 115 

I.F I MACHMX (J).LT.MCUT) GO TO 223 116 

IHTAGU) = 1 117 

GO TO 400 118 

223 R E P I J ) = 2. *PRBF*UAV*HCHAR/MU/R/ (T + 460. ) /2*681 1667 119 

IF < RHORF.NE.O. ) GO TO 221 120 
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IF ( ALPHA • NB. 0. ). AK1=-12.*MD0T< J) *MU*P1 / L/RHO 1/ 16 0. S3 33 3 121 

RHOREF = PRESS ( RDUF/2. )/R/ (T+460. 1*4.4756636 122 

GO TO 2 22 123 

221 RHOREF RHORF 124 

222 RER < J ) =i RHOREF*CAPV*HMEAN( J) /MU/12; 125 

Q(J) = 13'. 083**MD0T ( J ) 126 

K N ( J )« = 2.9 6* MA CHMX ( 0 ) /RE P ( J ) 127 

230 IF (TSKIP) GO TO 240 128 

POWER! J).= MU*AREA*CAPV**2/HMEAN< J)/6600. 129 

DELT J( J ) =.42 .42* POWER ( J ) /A8SI MDOT { J ) V/CP 130 

HTOTALtJJ = 42.42* POWER ( J ) 131 

TORUUE ( J)= POWER! J)*3. 3E4/ SPEED 132 

133 

DETERMINE FORCE AND CENTER OF PRESSURE 134 

135 

240 IF ( ALPHA. EQ.O. ) GO 10 250 136 

AK1 ~~i -L2.*M00T(J)* MU* PI / L/ RH01 / 160.83333 137 

K = 0 138 

KK = 0 139 

F(J)=i SIMPSUO. ,RDIF,PX,K)*L 140 

XC(J) = SLMPS1 (0**RDIF tPXXfKK) /F ( J)*L 141 

IF ( K.NE* 0) WRITE ' < 6 , 2 1) HME AN ( J ) 142 

IF (KK.NE.O) WRITE !6*24) HME AN t J )‘ 143 

GO TO 300 144 

250 F ( J ) = L*RD I F * ( 2 * *P L* < l.-C**3) /CC/3.-PREF) 145 

XC(J): = L*RDIF**2/F( J ) * ( 2 . *P 1 * ( . 4-C**3 ) /3 . / CC-PREF/2 * ) 146 

147 

DETERMINE PRESSURE, FILM THICKNESS, PRESSURE RATIO <P/P1>* 148 

AVERAGE VELOCITY, AND MACH NUMBER AT EACH GRID POINT 149 

150 

300 LF (ASKLP) GO TO 4C0 151 

P( J, 1) s PI 152 

IF ( ALPHA . L T. O'. ) P(J»L)=P2 153 

DO 320 1=1, NN 154 

IF ( UEQ.l) GO TO 310 155 

X ( J , L ) ^ FL0AT(L-1)*DELX+:X( J r l) 156 

H t J , I ). =i HME'ANl J)+ALPHA*(X( J*I )-RDIF/2. ) 157 

P t J, L) => PRFSS ( X! J , I ).) 158 

310 PRAT { J ^ I ) = P ( J , I )/PREF 159 

U A V R G ( J,.I )= ABS< MOOT! J ) ) *P1 /L/RH01/ *>( J *1 ) /HME AN ( J )/ 13. 405833 160 

MACH(JfI)= UA VR G ( J , I ) / A 161 

320 CONTINUE 162 

400 CONTINUE 163 

CALL STFN S S I HM E AN , F , I HT AG ,-S T I F F , N J ) 164 

165 

WRITE ROUTINE 166 

167 

600 WRITE 16,12) R , RHQ 1 ,.A , L , A RE A, SPEED , CAP V 168 

WRITE(6,.13) 169 

IF < /NOT. RSK'I P ) WRITE (6, 14) 170 

WRITE (6*55) 171 

IF (RSKI.P) GO TO 603 172 

WRITE (6*56) 173 

C 174 

603 DO 610 J=l, NJ 175 

WRITE (6,23) HMEAN(J) 176 

IF ( IHTAG! JV.E3.0) GG TO 609 177 

WRITE (6,20) 178 

GO TO 610 179 

609 WRITE (6, .15) MD OT ( J ) , Q ( J ); , M ACHMX ( J ) , REP ( J ) , RE R ( J ) * KN l J ) , F CJ > , 180 
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GO TO 610 


1 STIPF{ J) ,XC( J) 
IF (RSKIP) 

NORMALIZE F AND 


XC 


FBAR t J) =FU);/PDIJ=/RDIF/L 

XCBAR ( J ) =XC ( J) /RDIF 

WR I TE ( 6 ,.l 6 ) XCBAR( J) ,,F BAR ( J ) 

610 CONTINUE 

IF (TSKI.P) GO TO 612 

write (6,m: 

DO 611 J=1,MJ 

WRITE (6,23) HMEAN(J) 

I F ( LHTAGI J ) . EQ . O’) GO TO 613 
WR I T E ( 6 , 20 ) 

GO TO 611 

6 13 WRITE(6„22) POWER ( J )., H TOT AL ( J ) , OELT J ( J ) ».T ORQUE { J ) 

611 CONTINUE 

612 IF (ASKLP) GO TO 700 

WRITE X AND PRESSURE DISTRIBUTIONS 
JJ - 0 

DO 620 J = l, NJ 

IF ( LHTAGI J ) • ME . 0), GO TO 620 
JJ = JJ+1 

IF ( MOD (JJ,JMOD).EQ. 1) WRITE (6,57) 

WRITE (6,18) HMEAN(J) 

WRITE (6,19); ( X ( J , I ) ,.PRAT { J , I ) ,P( J, I) „UAVRG( J * I ) , MACH (I J, H, 
l I = UNN) 

620 CONTINUE 


PLOT ROUTINE 

700 KODE = 6 

P P ( 3 ). =. 0. 

P P { 4 ) = 0. 

DO 701 J ~ l , N J 

IF ( IHT AG ( J ) • NE . 0 ), HME AN < J ) --HM&AN ( J ) 

701 CONTINUE 

IF { ALPHA. EO.XT. ) GO TO 720 

CALL ARR.NG ( HME AN, XC *XPLOT , YPLOT , NJ ,NP ), 

PP( 1 1 NP 
WRITE (6,32): 

CALL PLOT XV tXP LOT, YPLOT, KODE ^PP) 

WRITE (6,40). 

710 IF (RSKIP) GO TO 720 

IF ( ALPHA. EQ.O. ) GO TO 720 

CALL ARRNGiHMEAN, F BAR, XPLOT, YPLOT ,NJ,NP): 
PP( 1 ). = NP 
WRITE (6,38) 

CALL PLOTXY (XPLOT, YPLOT, KODE,. PP> 

WRITE (6,421: 

720 IF (TSKIP) GO TO 730 

LF ( SPEED. EQ.O. ) GO TO 730 

CALL ARRNG ( HME AN, PO WER , XPLOT , Y PLO T , NJ , NP ) 

PP( 1 T = NP 

WRITE (6,34) 

CALL PLOTXY( XPLOT, YPLOT, KCDE,PP) 

WRITE (6,-43) 
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730 IF ( A SK I P ) GO TO 750 

DO 741 J*=1,NJ 

IF { LHTAGC J );.:NE.;0). GC TO 741 
WRITE (6*37): HMEAN(J) 

DO 740 1 = 1* NN 
XPT( L) =i X( J,.I.) 

YP T ( I ) ^ PRAT I J f I. ) 

740 CONTINUE 

CALL ARRNG < XP T * YP T * XP LOT YPL 0 T *.NN , NP ) 
PP( l ) => NP 

CALL PLOTXYIXPLOT* YPLOT , KCDE , P P ) 

WRITE (6,46) 

IF ( ALPHA • EQ ) GO TO 750 

741 CONTLNUE 

750 DO 751 J!=L*NJ 

IF ( IHT AG ( J ).'• NE*0 ) HME AN ( J ) =^-HME AN 1 J ): 

751 CONTINUE' 


CONVERT TO INTERNATIONAL UNITS AND PRINT 

800 I.F ( NOU I ) GO TO 100 
U I ( 2 1 ) s R*5* 3 8 0*9 5 
U I ( 2 2 ) = RHO 1*517. 2026 
UI (23) = A* ’• 3048 
UI (24)=. L*2.54E-2 
U I ( 2 5 ) =. ARE A*6 • 45 1 6 1— 4 
UI( 26) =. SPEED 1 / 1>0* 

UI (27)=, CAP V* • 3048 
WRITE (6*10) 

WRITE ( 6 *.54 V TITLE 

WRITE ( 6 * 5 0 )‘ (UI( I ) *1 = 1, 12) ,NN, (UI ( I) , 1 = 13* 15 ) , ASK I P*-RSK I P* 
1 TSKIP*,(UL( I.) *1=21 *27) 


WRI 

TE (6,13) 




IF 

< .NOT.RSKTP) 

WRITE 

(6 

, 14) 

WRI 

TE (6*51 ) 




IF 

( . NOKRSKJIP) 

WRITE 

16* 

56) 

DO 

810 J = 1 , N J 




UH 

= HMEAN(J)*2. 

54E-2 



WRI 

TE (6*23): UH 




IF 

( LHTAG ( J ) . E’Q* 

0) GO 

TO 

809 

WRI 

TE (6,20) 




GO 

TO 810 




809 UM 

= MDOT (J);** 755987E-2 



UQ 

= Q(J)**47194744 B— 3 



UF 

= F ( J ) *4 # 4482 2 1 6 



US 

= STLFF(J>*1. 

7512683E2 


UXC 

XC,( J)*2*54E-2 




WRITE (6*15) UM*UQ*MACHMX< J) *.REP< J ) * RER ( J ) * KN ( J ) , UF * US * UXC 
IF ( . NOT • RSK I P ) WRI TE<6 *16 ) XCBAR1 J) * F B A R ( J ) 

810 CONTLNUE 

IF (TSKIP) GO TO 812 
WRITE (6*52) 

DO 811 J = 1 * N J 

UH = HMEAN( J)*2.54E-2 

WRITE (6,23) UH 

IF( IHTAG(J) .NE-O) GO TO 813 

UP = POWER( J).*7.4569987E2 

UHT = HTGTAL(J)*i7. 597833 

UT = D E f L T J ( J ) / i • 8 

UTRK = TORQUE ( J 1*1 ♦ 3558179 
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WR I TE ( 6 »,22 ) UP,UHT,UT,UTRK 301 

GO TO 811 302 

813 WR I TE ( 6 ,.20) 303 

811 CONTINUE 304 

812 IF (ASKIPL GO TO 100 305 

JJ = 0 306 

DO 821 J=1,NJ 307 

LF { LHT AG { J V«*NE • 0 ) GO TO 821 308 

JJ = JJ+:i 309 

LF ( MOD ( J J » JMOD ) . E 0. 1 ) WRITE (6,57) 310 

UH = HMEANt J)*2.54E-2 311 

WRITE (6,53) UH 312 

DO 820 L-l,NN 313 

UX = X( J* L)*2.54E-2 314 

UP = P( J,I.)*6.8947572E3 315 

UV = UAVRGIJ, I ) *. 3048 316 

WRITE ( 6 ,.1 9 ) UXfPRAT<J,I). fUPfUV'tMACHi J , I ) 317 

820 CONTLNUE 318 

821 CONTINUE 319 

GO TO 100 320 

321 

1 FORMAT (13) 322 

2 FORMAT (6F12.6) 323 

3 FORMAT ( 1 2 A 6 ) 324 

10 FORMAT ( 1H1 ,.46HC0M PRES S LBLE SEALING DAM WITH SMALL TILT ANGLE) 325 

11 FORMAT ( lHOt-12-HlNPUT DATA - , / , 1H0 * 1 2HT I LT ANGLE - , F 8. 4 , 326 

1 9H RADIANS,,/, 1H0, 6X, 7HP2,PSI A,. 1 QX *7HP 1 # PS I A > 10X » 7HT * DEG F»10Xi 327 

2 2 OH V I SCO S I TY,LB-SEC/FT2, 10X, 16HM0LECULAR WEIGHT, 18X,5HCP/CVt/» 328 

3 1H ,.2G 1 7 • 5 ♦ F 1 3 • Of G25 • 5 » G3 1 • 5 ,028*5,/ , 1H0 4X , 9HR2 , INCHESt8X* 329 

4 9HR1,. INCHES, 9 X,8HL,. INCHES, 10X,15HRH0 t LB-SEC2/FT4,10X* 330 

5 20HRHQ(ROT) ,.L8-SEC2/FT4 , 10X, 17HN0 OF GRID POINTS,/, IH ,3G17w5, 331 

6 G21 .5,G27. 5,4.2 2, / , 1H0 , 6X , 5HN , RPM * 1 1 X , 8HV , FT/ SEC * 9X, 332 

7 l 5HC P , BT U/ LB- DEG R ,1 3 X , 6 H$K I P A,6X,6HSKIP R,6X,6HSK*P T*/,1H , 333 

R 2G17.5,,G18’. 5,L19,2L12) 334 

12 FORMAT { IHO, 17HBEGIN OUTPUT DAT A , / , IHO , 3X , 335 

1 30HGAS CONSTANT, FT-LB/LB (M)-DEG R,4X, 18HRH0( I) tLB-SEC2/F T4 i 6X , 336 

2 2 1H A ( SOUND SPE 60 ), , F T / SEC , / , 1 H , G24 ; 5 ,.G30 . 5 , G25 . 5 » / , IHO , 3X , 337 

3 8HL ^LNCHES , 14X ,.8H A?RE A ,IN2,14X, 9HSPFED ,*RP M , 13X t 8HV,FT/SEG*/, IH , 338 

4 G14.5t.3G22.51 339 

13 FORMAT ( 1H0^2X,9HMEAN F I L M , 4X , 6HM COOT ) , 6X , 1HQ *8X, 4HMACH , 5X , 340 

l 5HR E ( P )., 6X , 5HRE ( R ) ,.5 X , 7HKNUDSEN , 7 X , 1HF^7X,5HST IFF ,.7X t 2HXC ) 341 

14 FORMAT ( IH + *1 1‘OX, 2HXC , 8X *IHF ) 342 

15 FORMAT <lH+»llX v 2GlU3 v 2F9.3 v .?F10J3vFL2.3f2GlU3) 343 

16 FORMAT ( iH+,403X,2F10.3) 344 

17 FORMAT ( IHO, 16HMEAN F I LM LNCHE S ,.8X , 10HP0W ER , H . P . , 6X * 345 

1 18HSHEAR HFAT , BTU/MI N , 5X , 1 2HDEL ( T ) , DEG F* 8X , 12HT0RQUE , FT-LB) 346 

18 FORMAT ( IH0‘, l 1HMEAN FILM =,GI1.3,2X,6HINCHES, / , IHO, 5X^8 HX, INCHES,* 347 

1 4X*8HP/P<MIN) ,7X,5HP,PSI *6X,42HULAV) FT/SEC, 4X».7HMACH NO) 348 

19 FORMAT (IH ',614.3,4014.5). 349 

20 FORMAT ( 1H+ , 77X ,d 8HANALYS LS NOT VALID) 350 

21 FORMAT (HH0t26HF0R MEAN FILM THICKENESS =,G15.5,3X, 351 

1 14HF IS LNCORRECT) 352 

22 FORMAT (1H+, 15X,4G20.5) 353 

23 FORMAT ( IH , Gll-3) 354 

24 FORMAT (1H0,26HF0R MEAN FLLM THICKENESS =*G15.5,3X, 355 

I 15HXC I.S INCORRECT) 356 

30 FORMAT { 2HPT , 55HPL0T OF M(DOT) AND 0 VS HiMEAN) WHERE Q = 13J083*M 357 

UDOT) ) 358 

32 FORMAT ( 2BP T „2 3HP LOT OF X C C ) VS HIMEAN) ) 359 

33 FORMAT ( 2HP T„28HPLOT OF MOOT ( B AR ) VS HIMEAN) ) 360 



34 FORMAT ( 2HP T , 24HPLQT OF POWER V9 H ( MEAN ).). 361 

37 FORMAT ( 2HP T *.3'5HPL0T OF P/PIMIN); VS X FOR HIMEAN) =*G15.5> 362 

38 FORMAT (2HPT,.25HPLQT OF F ( BAR) VS HIMEAN).) 363 

40 FORMAT I2HPL..50HXC ON VERTICAL SCALE - HI MEAN ) ON HORIZONTAL SCALE 364 

1 * / « 2HPL4 /♦ 2HPL ,.10X , 56HXC IN INCHES - TO CONVERT TO METERS, MULTI'P 365 

2LY BY 2.54E-2 , / , 2HPL , / ,.2HPL , 10X ^61HH I ME AN ) IN INCHES - TO CONVER 366 

3T TO METERS', MULTIPLY BY 2.54E-2 ) 367 

41 FORMAT <2HPL,54HMI BAR) ON VERTICAL SCALE - HI MEAN) ON HORIZONTAL S 368 

1C AL E »/ ,.2'HPL', / » 2HPL * 10X,23HMIBAR)' IS 0 IMENS I ONLESS , / , 2HPL ,/ , 2HPL * 369 

2 10X ,6 1HHI MEAN) IN INCHES - TO CONVERT TO METERS* MULTIPLY BY 2.54 370 

3E-2) 371 

42 FORMAT 1 2HPL , 54HF I BAR ) CN VERTICAL SCALE - HIMEAN) ON HORIZONTAL S 372 

1C4LE,./».2HPL»/,2HPL» l OX ,23HFIBAR) IS DIMENSIONLESS, / *2HPL,'/» 2HPL » 373 

2 LOX *6 1HH ( MEAN I IN INCHES - TO CONVERT TO METERS* MULTIPLY BY 2.54 374 

3E-2) 375 

43 FORMAT 1 2HPL ,.5 3HP0WER ON VERTICAL SCALE - HIMEAN! ON HORIZONTAL SC 376 

1ALE, /,2HPL, /..2HPL, 10X.61HP0WER IN HORSE POWER - TO CONVERT TO WATT 377 

2S , MULTI.PL Y BY 745 . 7 , / , 2HPL > / , 2HPL , 10X , 6 1HH I MEAN ) IN INCHES - TO C 378 

30NVERT TO METERS,. MULTIPLY BY 2 . 54E-2 ,./ , 2HPL , / , 2HPL , 50HFOR SHEAR H 379 

4E AT IN BTU/MIN, MULTIPLY POWER BY 42.42) 380 

46 FORMAT ( 2HP L „50HP/P I M I.N ) ON VERTICAL SCALE - X ON HORIZONTAL SCALE 381 

1> / * 2HPL *./■, 2HPL , 10X ,25HP/P(MIN) IS D I MENS I ONLESS ,/»2HPL»/,2HPL»10X* 382 

2 55HX IN INCHES - TO CONVERT TO METERS*. MULTIPLY BY 2.546-2) 383 

50 FORMAT UH0,12HTILT ANGLE =,F8.4*9H RADI. ANS, / , 1H0, 5X , 7HP 2, N/M2, 384 

1 9X, 7HPl*N/M2, 10X,7HT,0EG K»10X*.16HVISC0SITY,N-S/M2*1QX» 385 

2 16HM0LECULAR WEIGHT, 1 4X ».5HCP/CV*./ , 1H , G14 . 5 , Gl 7. 5 ,.F 14. 0,G24. 5, 386 

3 G26.5,G2?.5,/,lH0,4X,9HR2,METERSn7X,.9HRl , METERS, 9X,8HL, METERS* 387 

4 13X^9HRH0,K'G/M3,12X.,14HRH0IR0T):,KG/M3,iaX, 388 

5 17HN0 OF GRID POINTS, /,1H , G 14. 5 *Gl 7.5 „G 18 . 5 ,G21 . 5 , G22 . 5 , 1.21 , / * 389 

6 IH0,6X,.5HN,.RPS,12X,5HV,M/S,9X,13HCP,J/KG-DEG K, 390 

f 16X,6HSKjIP A ,.6X , 6H SKI P R,6X,6HSKIP T,/,.IH , G15 .5 * G 17 . 5, 391 

8 G18.5,.L?0,2L12,./,1HO,17HBEGIN OUTPUT OAT A, / , 1H0* 3X* 392 

9 23H0AS CONSTANT, J/KG-DEG K,10X,.12HRH0( 1) ,KG/M3».10X* 393 

X 1 8H A ( SOUND SPEED )., M/S ,/, 1H , G20. 5 , G3 1 . 5 , G2 3. 5 , / , IHO* 3X, 394 

1 8HL».METERS*.15X»7HAREA,M2».15X ».9HSPEF0 , RPS , 1 5X ,5HV *M/9</*lH *G14.5» 395 

2 3G23.5), 396 

51 FORMAT UH ».4X*6H METERS, 5X»6HK,G/S£C*5X».4HSCMS ,6 X*5H (MAX) ,25X* 397 

1 6HNUMBER.5X ,.7HNE W TONS ,4X,4HKG/M,6X,6HMETERS) 398 

52 FORMAT I1H0.16HMEAN F I LM ,ME TERS , 7X, 1 1HP0WER ,.W AT TS * 8X , 399 

1 16HT0TAL HEAT , WATTS *-6 X, 12HDEL IT) , DEG K».8X, 10HT0RQUE ,'N— M > 400 

53 FORMAT I IHO „1 1HMEAN FILM = ,G 14. 5 * 2X ,6HMET ERS, / , IHO , 5X, 8HX , METER S,i 401 

1 4X, 8HP/P(MI.N).,7X ».6HP ».N/M2*6X, 11HUIAV). *.M/SEC, 5X,7HMACH NO) 402 

54 FORMAT UH0..12A6) 403 

55 FORMAT I 1H , 3X , 6H I NCHE S*6X ,6HLB/MI N, 5X ,-4HSCFM ,6X , 5H I M AX ) * 25X* 404 

1 6HNLMBER ,8X».2HLB»6X»5HLB/IN»5X»6HINCHES) 405 

56 FORMAT I IH + » 109X, 3HB>AR»7X»3HBAR) 406 

57 FURMAT (1H1) 407 

C 408 

END 409 
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$«BFTC PRESR 
C 

C PRESSURE FUNCTION 

C 

FUNCTION PRESS < X) 

COMMON/ I N TSRJL/PLt P2 , ALPHA ^C^RDIF ,HM,AK 
C 

YY = X 

LF ( ALPHA*NE.O. ) GO TC 100 

Y = P1*SQRT{ 1.-C*YY/RDIF) 

GO TO 101 

100 AA = HM-ALPHA*RDI F/2. 

YD I F = YY/ROIF 

q = Pl**2*( l.-YDIF*C)+AK/ALPHA*< ( 1 . -YD I F ), /A A* #2-< l./C 
1 ALPHA*YY+AA ) **2-YDLF/ < ALPHA*RD I F + AA ) **2 ) ) 

Y « SQRT(Q) 

101 PRESS =n Y 
RETURN 
END 


410 

411 

412 

413 

414 

415 

416 

417 

418 

419 

420 

421 

422 

423 

424 

425 

426 

427 

428 


SI.BFTC FPDX 429 

C 430 

C EXTERNAL FUNCTION FOR INTEGRAL tP~Pl) DX 431 

C 432 

FUNCTION PX<Y). 433 

COMMON/ 1NTGRL/P1 » PREF * A ,B *0 t D t*E 434 

YY=Y 435 

PX = PRESS ( YY) -PREF 436 

RETURN 437 

END 438 


$ 1 B F T C FPXDX 

439 

C 

440 

C EXTERNAL FUNCTION FOR INTEGRAL IP-PDXDX 

441 

C 

442 

FUNCTION PXX(Y) 

443 

COMMON/ I.NTGRL/P1»PREF ? A f B*C,D,E 

444 

Y Y = Y 

445 

PXX ^ t PR'ESSt YY)-PREE) *YY 

446 

RETURN 

447 

END 

448 
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IBFTC DERLV 


LAGRANGE* NUMERICAL 


DIFFERENTIATION OVER MAXIMUM OF 5 


POINTS 


SUBROUTINE STFNSS < XXVYY».ITAG*DDYfrMAX) 

DIMENSION XX(50)»YY<50) 9 ITAG<50);tDDY( 50) f X(50l;f Y(50! f DY(50) *A(5t9> 


ELIMINATE I NVA'L 1.0 POINTS AND ARRANGE VALID POINTS I.N ASCENDING 
ORDER. IF THERE ARE LESS THAN 2 VALID POINTS* NO 0 I F FERENT I AT I ON 
LS POSSIBLE; 


LOO MM = 0 

DO 110 It MAX 

IF ( I TAG ( M ) . NE 1 . 0 ) GO TO 110 
MM S MM+:i 
X(MM); = XX ( M ) 

Y ( MM ); = Y Y ( M ); 

110 CONTINUE 

LF (MM.GE.2) GO TO 130 
DO 120 M= 1 , MAX 
120 DDY(M) =. 0. 

RETURN 

130 CALL SORTXYlXt Y»MM) 

SET UP MATRIX OF X DIFFERENCES FOR EACH POINT X ( K ) 

200 N = M LNOi MM f ,5 ) 

DEBUG N f MM 

DEBUG <X(M), Y ( M) * M= 1 , MM ) 

DO 250 K=1,MM 

1ST == MAXO ( Kl— 2 , 1 L 

I.S T MLNOIMM-N+l, 1ST) 

IN = IST+N-1 

DEBUG K!f.LST f IN 

DO 2 LI LI.- 1 ST f IN 

I = LI-IST+1 

DO 2 LO J J’= I ST t I N 

J = JJ-IST+1 

A( I f J) = XIII) -X( JJ‘) 

2 LO CONTLNUE 

DEBUG ( A ( L, J J ) t JJ=1,N) 

211 CONTLNUE 

FORM SUMS AND PRODUCTS FOR DERIVATIVE FORMULA 


220 


230 

231 


SI 

= 0 *. 




S2 

= 0 . 




P 2 

* 1. 




DO 

231 LI - 1 

ST, I 

N 


IF 

< IL.EQ.K 

) 

GO 

TO 231 

I = 

LL-LST + 

l 



PI 

= X< I u- 

X { K ) 



S2 

= S2-1./ 

PI 



P 2 

= P2*Pl 




DO 

230 j:j = l 

t N 



IF 

I l.NE; J J 

); 

Pl = 

=Pl*A(l # JJ) 

CON 

ITLNUE 




DEBUG Pl,Yl 

ii ) 



SI 

^ Sl+YU 

D;/p 

1 


CONTINUE 
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470 
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480 
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4 85 

486 
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488 
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IF I (N/2)*2.NE.N) S2=-S2 509 

DEBUG S 1 V 52 ».P 2 510 

511 

DERIVATIVE 512 

513 

KY = S2*Y(K)+:P2*Sl 514 

DY(K). =. KY 515 

250 CONTLNUE 516 

517 

PUT CALCULATED DERIVATIVES IN ORDER TO CORRESPOND TO INPUT XX 518 

ARRAY- 519 

520 

300 DO 320 M=. 1 » M AX 521 

IF 1 LTAG1M) LNE.O) GO TO 320 522 

DO 310 1.1 = 1, MM 523 

I.F I XXI M) .NF. XI II.) ) GO TO 310 524 

IF I (N/2),*2.N6.N) GO TO 311 525 

DDYIM) = -D Y ( I.I ) 526 

GO TO 320 527 

3 1 1 DDY ( M) =5 DY1I I ) 528 

GO TO 320 529 

310 CONTINUE 530 

320 CONTLNUE 531 

DEBUG (XIM) ,DY(M) ,M=1,MM) 532 

533 

RETURN 534 

END 535 


TC SMPSRl 536 

537 

NUMERICAL INTEGRATION BY SIMPSONS RULE 538 

539 

FUNCTION SI MPS l (XMIN*XMAX »FUNC 1 t-KER )' 540 

DIMENSION VI 200) ,.H ( 200 ) , A ( 200) ».B { 200 ) ».C 1 200 ) » PT 200 ) v E 1200 ) »NE ( 200 ) 541 

EQUIVALENCE l E , NB) , < TEST „NTEST ) 542 

T=3 . OE- 5 543 

V11MXMIN 544 

11(1) =>0. 5*( XMAX-XMl N ) 545 

All) =FUNC L ( XM I.N ) 546 

B< IMFUNCKXMIN+Hdl ) 547 

C( l ) =;FUNC1( XMAX) 548 

PI 1 ) =.H ( 1)*(A(1)+'4,0*B( l ) +'C ( 1 ).) 549 

E(1)=5P(1) 550 

AN S = P ( 1 ) 551 

N= 1 552 

FR AC =2 . 0* T 553 

1 FRAC=sO. 5*FR AC 554 

2 TEST=;ABS(FRAC*ANS) 555 

K=N 556 

3 DO 7 L=.1,K 557 

4 IF (NTEST-IABSINEI I)) ) 5,5,7 558 

5 N = N+l 559 

V ( N ) sV ( I ) +H < I ) 560 

HIN)40.5*H( I). 561 

A 1 N ) =;B ( I .) 562 
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BIN) ^FUNC L ( V ( N ) +:H ( N ) ) 563 

C(N)=iC(n 56 A 

P I N ) ^H< N) * ( A ( N ) +4* 0*B( N ) -*iC ( N). J 565 

Q=P(I) 566 

H(I)=5H(N) 567 

BII) =iFUNCl( VI I.)+:H( I ).) 568 

C I I ) =sA I N) 569 

pm =;H( u*(A(u+:4iO*B( n+:c i m 570 

Q=P( L)+:P(N)-Q 571 

ANS= ANS+Q 572 

ElIhQ 573 

EtN)=iQ 574 

6 IF IN-200) 7*13,13 575 

7 CONTINUE 576 

8 IF IN-K) 9,9*2 577 

9 Q = 0.0 578 

10 DO 11 1^1, N 579 

11 Q=Q*E< I ) 580 

12 IF (ABS(U)-T*ABS14NS) ) 14,14,1 581 

13 KER = KER+:i 582 

14 AN S=0. 0 583 

15 DO 16 I=il,N 584 

16 ANS = ANS + P I I > 585 

SIMPS1=>IANS+'Q/'30 I . 0) /3 . 0 586 

17 RETURN 587 

END 588 


SIBFTC RRANG 589 

C 590 

C SUBROUTINE TO ARRANGE ARRAYS TO BE PLOTTED 591 

C 592 

SUBROUTINE ARRNG I X , Y , XP , YP , N , I): 593 

DIMENSION XI25) ,YI25> *XPI25) ,YP(25) 594 

I. = 1 595 

DO 100 J= 1 , N 596 

IF I XI J ULT.O. ) GO TO ICO 597 

X P I I ). => Yl JV 598 

YP I I 1 = XI J): 599 

I = L+l 600 

100 CONTLNUE 601 

I = 1-1 602 

CALL SORTXY (XP, YP,I.) 603 

I.I = I./2 604 

DO 101 J*=l,l.I 605 

IJ = L-J+l" 606 

T = YPIJ) 607 

YPI J ); =. YP ( I.J) 608 

YP I I J) =; T 609 

T = XPIJ) 610 

XPIJI s XPIIJI 611 

XPI I J) = T 612 

101 CONTLNUE 613 

RETURN 614 

END 615 

>01* UNIT05, EOF. 
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APPENDIX D 


EXAMPLE PROBLEM 

An example of the use of the computer program is given with the following condi- 
tions: Air, internally pressurized at 65 psia, is to be sealed from ambient air at 15 psia. 
The mean temperature is 100° F. The sealing dam outside diameter is 6. 630 inches, 
inside diameter is 6. 530 inches, and the design speed is 1398 rpm. Mean film thick- 
nesses of 0. 1 to 1. 0 mil are to be investigated in increments of 0. 1 mil for parallel 
surfaces and a positive and negative tilt angle of 1. 0 milliradian. Thus, 


Pressure at inner radius or inlet, P^, psia 

65 

Pressure at outer radius or outlet, P 2 , psia 

15 

Temperature, T, °F 

100 

Inner radius, Rp in. 

3. 265 

Outer radius, R 2 , in. 

3.315 

Mean film thickness in increments of 0. 1 mil, h m , mils 

0. 1 to 1. 0 

Relative inclination angle of sealing dam surfaces, rad 

0, ±1.0xl0" 3 

Design speed, rpm 

1398 

Molecular weight of gas, m, lbm/lb-mole 

28.966 

Absolute or dynamic viscosity, \i , lb -sec/ft 2 

3.96X10" 7 

Specific heat at constant pressure, C p) Btu/(lb)(°R) 

0. 24 

Specific heat ratio, C p /C v 

1.4 


The data sheet for this sample problem is shown in table in. The sample output with 
all possible output options is shown in both English units and International System of 
Units. Note that the analysis is not valid for a mean film thickness of 0. 4 mil and larger. 
The Mach number has exceeded 1/ y 'y = 0. 845. The execution time for this problem is 
about 0. 39 minute on the Lewis computer. 
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05 


LUMP REbb i dL t bt A L I iMG DAM irflTH SMALL JILT ANoLE 


bAMPLb PKUbLcM 
INPUT L>a T A - 


TiLT ANbLt - J. KAJIrtNS 


P2,Pbi A 


Pl,PbU 

T » Ob G P 

VISCUSITY, LB-SEC/ FT2 

MOLECULAR 

HEIGHT 


CP/ CV 

IS.OOoC 


63. UU 00 


100 . 

0. 

39600 E-06 

28.9660 


1.4DDD0 

R2» INLHE b 

Ki 

, INUHLb 

L » I NChb b 

RHO.LB 

-SEC2/FT4 

RHO(ROT) ,Lb-SEC2/FT4 

NO 

OF GRID POINTS 

3.3i5uC 


3.26300 

20. 

0000 

0. 

93000 E-02 

0.600006 

-02 


11 

N , KPM 

V 

r , P T/ btU 

CP,8TU/LB-JEG K 


SKIP A 

SKIP R SKIP 

T 



U 


4C. 0000 

0 . 

24000 


T 

F F 




BfcC-lN OUTPUT 

OATA 










GAb UUNblANl ,FI-Ltt/L8lM)-0Eb K 

RHUUJ 

, L6-SEC2/ FT4 

Al SOUND SPEED), FT/SEC 





33 .3322 



0 . 93000 £—0 2 


1160.08 





L v iNLHEb 


AREA * I N2 


b PEED , RPM 

V, FT/SEC 





Lu .GOOD 


1.0333d 


1393.21 


40.0000 





M t AN MLM 

M 10UF ) 

0 

MACH 

Kb(P) 

Rh(R) 

KNUDSEN 

F STIFF 

XC 

XC 

F 

INUHLb 

Lo/MiN 

SUP M 

(MAX) 



NJMBER 

LB LB/IN 

INCHES 

BAR 

BAR 

C. luOE— 02 







ANALYSIS NOT VALID 




C.9U0L-03 







ANALYSIS NOT VALID 




C. bG0t-03 







ANALYSIS NOT VALID 




C.7CUL-03 







ANALYSIS NOT VALID 




C.600E-03 







ANALYSIS NOT VALID 




C • 500b- U 3 







ANALYSIS NUT VALID 




C.45UE-G3 







ANALYSIS NOT VALID 




C.4G0E-03 







ANALYSIS NOT VALID 




C. 35GE- 03 

0.1 Oo 

2. 1 7*t 

0.711 

273.048 17.677 

0.008 

30.208 -0 

0. 170E-01 

0.339 

0.604 

C.30ub-U3 

0. X 03 

1.369 

0.522 

171.949 15.152 

0.009 

30.208 0 

0. 170E-01 

0.339 

0.604 

C.23UE-U3 

u. 6U5E-G1 

0. 792 

0.363 

99.307 12.626 

0.011 

30.208 0 

0. 170E-01 

0.339 

0.604 

C • 2U0E- U 3 

0.3 lot- Ci 

C.40o 

0.232 

50.948 10.101 

0.013 

30.208 0 

0. 170E-01 

0.339 

0.604 

C* 130E- u3 

0. 13 lb- 01 

0.1 71 

0.131 

21.494 

7.576 

0.018 

30.208 0 

0. 170E-01 

0.339 

0.604 

C. 1U0E-03 

0. 387t- C2 

0. 3 0 7b -01 

0.058 

6. 36 8 

5.051 

0.027 

30.208 0 

0. 170E-01 

0.339 

0.504 

MEAN FILM , AN Uric S 

POwt R,H .P. 

SHEAR HEAT ,BTU/MIN 

DEL ( T ) , DEG F 

TORQUE, FT-LB 




C. 100L-02 







ANALYSIS NOT VALID 




C.900E-03 







ANALYSIS NOT VALID 




C.80uE-03 







ANALYSIS NOT VALID 




C. 700E-03 







ANALYSIS NUT VALID 




C.6GUE-03 







ANALYSIS NOT VALID 




C.500t-U3 







ANALYSIS NOT VALID 




C.450E-03 







ANALYSIS NOT VALID 




C.40UE-03 







ANALYSIS NOT VALID 




C.350E-03 


0.26330E-03 


0.12026E-01 


0.30161 

0. 67150E— 02 




C.300E-U3 


0. 33 073b— 03 


0. 14030E— 0 1 


0.55877 

0.78342E-02 




G.250E-U3 


G. 396.906—03 


0.16836E-01 


1.15867 

0. 940 I0£— 02 




C.200E-U3 


G.49612E-03 


0.21 045 E— 0 1 


2.82878 

0.1175 IE— 01 




C. 150E-U3 


0. 66 1 49E-G3 


0.2806 IE-01 


8.94033 

0. 15668E — 01 




C.100E-03 


0.99224E-U3 


0.4209 IE-0 1 


45.2604 

0.23 503E— 01 





PLOT OF POWER V 5 HIMtANJ 


U. 000 100 

0, 001 00* 

1 

1 

1 

1 

1 

1 

1 

1 

1 

O.OOOdUl 

1 

1 

1 

1 

1 

1 

l 

1 

1 

U. QUGbOi 
1 
1 
1 
1 
1 
1 
1 
1 
1 

C. 000401 
1 
1 
l 
1 
1 
i 
1 
1 
1 

C. 000201 

u.OOOl uo 


0.0u0125 0.000150 0.300175 0.000200 0.000225 0.000250 0.000275 0.000300 0.000325 


0.000125 0.000150 0.J0U175 0.000200 0.000225 0 .000250 0.000275 0.000300 0.000325 

PUwfcR UN VERTICAL 5CAl c - HlflEAN) UN HURlZuNfAL 5 C ALE 

PCWEK IN HUKSL PUwER - Tu CONVERT TU WAT 1 5 f MULTIPLY BY 745.7 
HMtANJ IN INCHED - fU LONVEKT TU MtTERS, MULTIPLY BY 2.54E-2 


0.003350 


* 


0.000350 


FUK 5HEAK HEAT IN b TU/Ml N» MULTIPLY POWER BY 42.42 


CGrtPKEbalBLt SEALING DAM wllrt aMAi . l TILT ANGLE 


SAMPLE PRGtJLtM 

TILT ANGLE = O. RADlrtNS 


P 2 , n /M 2 

PI ,i\t/M2 

T» JcG 

K 

VlaCGSITY, N— S / M 2 

MOLECULAR WEIGHT 

CP/C V 


U .10342t 

Cto U • 44 bl ot ub 

311 

* 

0. 1896 1 E— 0 4 

28.9660 


1.40000 


K^frthl tKi Kl,MtTERS 

L, METERS 

KHOi 

i K G/M3 

RHO ( ROT ) , KG/M3 

NO 

OF GRID POINTS 

0 .8420 lt~ 

Cl 0.82 93 It — 01 

U.508O0 

4.80998 

3. 10322 


11 


N,KPa 

V,M/S 

CP, J/KG 

-OfcG K 


SKIP A 

SKIP R SKIP T 




0 

12. 1920 

1004 . 

78 


T 

F F 




BEGIN JUlPuT 

OA TA 









GAS LU.\lalANT,J/KG-UtG K 

KHC( 1) ,KG/M3 

AtSGUND SPEED), M/S 




2d7 

• Gob 

4. 80 99 d 


353.591 





L »ME 1 tKa 

AKtA ,M2 


SPEED 

, KPS 

V,M/S 




u .soaoo 

0* 66b8aE~ 

•03 

28.2202 

12. 

1920 




MtAN FILM 

M( 0 U T i 0 

MALH 

Kt(P) 

REIR) 

KNUDSEN 

F STIFF 

xc 

XC 

F 

METEKa 

KG/ atC aCMa 

(MAX) 



NUMBER 

NEWTONS KG/M 

ME TERS 

BAR 

BAR 

G.254E-U4 






ANALYSIS NOT VALIO 




0.229E-U4 






ANALYSIS NOT VALIO 




U. 203E-04 






ANALYSIS NOT VALIO 




0. i 78 1— 04 






ANALYSIS NOT VALIO 




0.152E-U4 






ANALYSIS NOT VALID 




0.127E-U4 






ANALYSIS NOT VALID 




G. U4E-04 






ANALYSIS NOT VALID 




0. 102E-U4 






ANALYSIS NOT VALIO 




C.889E-U5 

O.i2oc-G2 0. lu3E-02 

0.7 11 

273.048 

l 7.677 

0.008 

134.373 -0 

0.43LE-03 

0.339 

0.634 

G.762E-05 

O.791E-03 O.b46t-03 

0.322 

171.949 

15.152 

0.009 

134.373 0 

0.431E-03 

0.339 

0.604 

0. t35fc~0i> 

0. 458E- 03 0.374t-03 

0.363 

99.507 

12.626 

0.011 

134.373 0 

0.431E-03 

0.339 

0.604 

u • 5 cat - ud 

0.2 34t— 03 0.191E-03 

0.232 

50.948 

10.101 

0.013 

134.373 0 

0. 431E- 03 

0.339 

0.604 

C • 38 IE- 05 

U . 9 89E— 04 0.o07fc-04 

O.lii 

21.494 

7.576 

0.018 

134.373 0 

0.431E-03 

0.339 

0.604 

0 • 254E- 05 

0.293E-04 0.239E-04 

0 .0a8 

6.36 8 

5.051 

0.027 

134.373 0 

0. 431E-G3 

0.339 

0.604 

MEAN F ILM , Me 1 LR a PUwEK, wATTS 

TuTAL HEAT, 

WATTS 

D£L( I ) , DEG 

K TORQUE ,N-M 




U.^54E-04 






ANALYSIS NOT VALID 




0.229E-04 






ANALYSIS NOT VALID 




0 • * Q3E- 04 






ANALYSIS NOT VALID 




C • 1 7at- 04 






ANALYSIS NOT VALID 




0. 152E-C4 






ANALYSIS NOT VALID 




0. 127E-04 






ANALYSIS NOT VALID 




0. U4E-04 






ANALYSIS NOT VALID 




G.1G2E-04 






ANALYSIS NOT VALID 




t.869E-Q5 

0.21140 


0.21163 


0. 16756 

C.91043E-02 




0.762E-05 

0. 24ob4 


0.24690 


0.31043 

0 . 10 622E-0I 




C.635E-G5 

0.29397 


0.29628 


0.64370 

0.12746E-0I 




0 • S'CbE- 05 

0. 36996 


0.37035 


1.57154 

0.15933E-01 




0 • 38 IE- 03 

0.49328 


0.49381 


4.96 685 

0. 21243E— 01 




0 • 254E- 03 

0. 73 991 


0.74071 


25.1447 

0.3186 5E-01 








CUMPKtSSIbL t BALING DAM wi TH 8 MALL TILT ANuLt 
SAMPLt PROBLEM 
INPUT DATA - 

TALI AN Out = J . RAUIAnS 


P 2, PS I A P 1 » P Si A T » DtG F V IS COS l T Y ,L 8-SEC/ FT 2 MOLECULAR WEIGHT CP/CV 

Id.UuUL 08.UOUO UO. G.396G0E-06 28.9660 1.40000 


K2, iNCHtS 

. K1 

, 1 Nutifc S 

L, INCHES KHU,L8 

-SEC2/FT4 

RHO(ROT) ,LB-S£C2/FT4 

NO 

OF GRID POINTS 

3 • 3 i5U C 

3.2o3UU 

U 

0 


0 



II 

N,rtPM 

V 

,P i/sfcb 

CP,bl U/Lts-OtG R 

SKIP A 

SKIP k SKIP 

T 



1398, uC 

0 

U. 

24000 

F 

F F 




BEGIN OUTPUT 

OAT A 









GAS CONS TAN T ,F 1-LB/Lb i M) -ufcb K 

KHUIli 

, LB-S LC2/ FT 4 AtSOUNO SPEEO), 

FT/SEC 





33.3322 



U . 97 3 7 i h-02 

1160 .08 





L , INCHts 


AREA, IN2 


SPEEO, RPM 

V, FT/SEC 




20.671/ 


1. 03336 


1393.00 

40. 1375 





Mt AN FILM 

MIDOT i 

W 

MACH 

KE(P) HEIR) 

KNUOStN 

F STIFF 

XC 

XC 

F 

INCH to 

LO/MIN 

SoFM 

i MAX) 


NUMBER 

LB LB/IN 

INCHES 

BAR 

BAR 

C. 10UE-02 






ANALYSIS NOT VALID 




C.9UOE-03 






ANALYSIS NOT VALID 




C.ti0Ufc-03 






ANALYSIS NOT VALID 




C • 7U0E-U3 






ANALYSIS NOT VALID 




C.600E-U3 






ANALYSIS NUT VALID 




C.5CUt-U3 






ANALYSIS NOT VALIQ 




C.45uE-03 






ANALYSIS NOT VALID 




C. AGUE- Q3 






ANALYSIS NOT VALID 




G.j50E-u3 

u. 180 

2. 332 

0./ U 

273.048 20.869 

0.008 

31.223 -0 

0.I70E-01 

0.339 

0.604 

C.3CUE-03 

U. lii 

1.4dl 

u .322 

171.949 17.905 

0.009 

31.223 0 

0.I70E-01 

0.339 

0.604 

C. 250E- U3 

U. c 33c- Cl 

U. 7 

U. 363 

99.5 0? 14.921 

0. Oil 

31.223 0 

0. 170E-01 

0.339 

0.604 

C.20uE-uJ 

U.333E-U1 

0.439 

0.2 j 2 

30.94b 1L.937 

0.013 

31.223 0 

0. L70E-01 

0.339 

Q. 534 

C.15UE-G3 

0. 142E-G1 

U. 133 

0.131 

21.494 8.953 

0.018 

31.223 0 

0.170E-01 

0.339 

0.504 

G. IOUfc-03 

0.4 19E- 02 

U.549E-01 

U.OSo 

6.368 5.968 

0.027 

31.223 -Q 

0. 170E-01 

0.339 

0.604 

MEAN FILM, INCHES 

PUrtfc K ,H . P. 

SHEAR HEAT , BTU/M IN 

DEL { T ) , OEG F 

TORQUE, FT-LB 




G. 10UE-02 






ANALYSIS NOT VALID 




C.9U0E-U3 






ANALYSIS NOT VALID 




C • ttOUfc— U 3 






ANALYSIS NOT VALID 




C.7GOE-U3 






ANALYSIS NOT VALID 




C.6CUE-03 






ANALYSIS NOT VALIO 




G.50GE-03 






ANALYSIS NOT VALID 




C.450E-U3 






ANALYSIS NOT VALID 




C.40UE-U3 






ANALYSIS NOT VALID 




C.35UE-03 


0.28S45E-Q3 


0.12109E-01 

0.28063 

0. 67381E— 02 




G.30UE-U3 


0. J3302E-Q3 


0 . 14 12 7 b-0 1 

0.51990 

0.7861 IE-02 




C.25UC-03 


0. 39 9o3E -U3 


0.169P2E-U1 

1.07807 

0.943336-02 




C. 200E- 03 


0.49954E-03 


0.21 190 1-0 1 

2.63200 

0.11 792E-01 




C.150E-U3 


0 . 6o603E— 03 


0.28234E-01 

8.31843 

0 • 15722E—Q1 




C.100fc-U3 


0.99907E-03 


0 .4238 1E—0 1 

42.1120 

0. 23583E-01 





CO 


MEAN FILM 


0.3;>UE-C3 INCHES 


X» INCHES 

P/P1MI Ni 

P,PSl 

U (AN/ > ,FT/StC 

MACH NU 

0 

4. 3333;> 

o5. OvJUO 

190 .363 

0.16410 

C .5 UOd- C2 

4.12311 

ol. 64 oo 

20U .07 1 

0 . 17246 

G. lo je- C L 

3.9015/ 

36. 3233 

211.432 

U .18226 

0 • i30c- L l 

3 • 6060 7 

3 3. 00 UU 

224.977 

0.19393 

U.2u0t- Cl 

3 . 4 1 :>u5 

31.2348 

241.511 

0.208 18 

C.250E- Cl 

3. I44c6 

47. lo 99 

262.323 

0.226 13 

0.3uUc- C i 

2 .04600 

42. 72 00 

2 d 9 • 6 4 / 

0 . 24968 

C.350E- Cl 

2 • j lool 

37. 7492 

327.788 

0.28256 

C.400c- Ci 

2.13437 

32.0136 

3do .4 90 

0.333 16 

u **t30c- c 1 

1 .0660 7 

23.0000 

494 .949 

0.42665 

0. 6U0E- C 1 

1.O0OUO 

13. OuOO 

824.915 

0.71109 


MtAN Pil_M = J.jU0t-03 INCHES 


X, INCHES 

P/PlrtiN) 

P, PSi 

UiAVJ .FT/dEC 

MACH NU 

0 

4. 33333 

65. 0000 

139.860 

0 . 120 56 

0. 50ut- C2 

4.12311 

ol. 64 06 

146.991 

0.126 71 

U.iOOE-Ci 

3.9C157 

56. 5235 

133.338 

0.13390 

0. 130E- Ci 

3 . 6o667 

33.0000 

165.289 

0.14248 

0.2UOE- Cl 

3 .4 1 3o5 

51. 2346 

177.43o 

0.15295 

0.230E- Cl 

8 • 14466 

47.1699 

192.727 

0. 16613 

0 • 300c“ C 1 

2 . 84du0 

42. 7200 

212.802 

0 . 18344 

0.350E- Cl 

2.51661 

3 7.74 92 

240.824 

0 .20739 

U.4u0c- C 1 

2.13437 

32. 01 5o 

283.932 

0.24477 

0.430E- Cl 

1.66667 

25. OuOO 

363 .636 

0 . 3 1 3 46 

U .5cOE- C 1 

1.00000 

13. 0000 

6 06 .OoD 

0.52243 


MEA\i FIlM = O.250E-03 INCHES 


x # Inches 

P/P tMINJ 

P,PS1 

J(AV) , FT/SEC 

MACH NO 

u 

4.33333 

6 3.00 UO 

97.1251 

0 *83723 E— 01 

0.50UE-C2 

4.UJU 

61 . 84 66 

102.077 

0. 87992 E- 01 

u. lOOt-Cl 

3.9C157 

58.5235 

107.8 73 

0.9 29 88 E- 01 

0.15uE-Ci 

3 . 6666 7 

55. 0000 

114.784 

0.989 45 E— 01 

U.200E- Ci 

3.41565 

51.2348 

123.223 

0.10622 

C .230e- L 1 

3.14466 

4 7. 1 o 99 

133.838 

0.11537 

C.3U0E- C l 

2 .648U0 

42. 7200 

147.779 

0.127 39 

0.35UE- Cl 

2.51661 

37. 7492 

167.239 

0. 144 16 

0 • 4uOE- Ci 

2.13437 

32.0156 

197.189 

0.16998 

0.430E- C 1 

1 . o6b6 7 

25.0000 

252.325 

0.21768 

0.5 uUE“ C 1 

l.uOOOO 

15. 0000 

420 .875 

0.36280 



Mb Aft FILM 


0.200E-C3 INbHti 


Mb AN 


MtAft 


A i 1 l>ib He 8 

P/P (MiftJ 

P i P 3 I 

0 I A V ) »F T/ 3 bb 

MACH NO 

0 

4. 83.333 

63. 00 Ob 

ol . 1600 

0 .53583E-01 

l.3Uub - it 

4-. 12311 

6 1 . 84 D 6 

65.3294 

0 .563 15E-OI 

0 a 1 jUt- Cl 

3 . 9013 7 

38. 5235 

69.0390 

0 .595 12E-01 

0. i50t- Ci 

3 • 0660 7 

35 . 0000 

73 . 46 19 

0 . 63325 E-O t i 

U. 200 c- Ci 

3.4i5.>5 

51.2348 

78. 86 0 o 

0 .679 79 E - 01 

b.^3Ub- C 1 

3 • i44o6 

4 7. 1099 

83.8364 

0 . 73837E-01 

0 . jUuc- C i 

2 • 04 abb 

42. 7200 

94.3/8/ 

0 .815286-01 

G.33Ut- Cl 

2.5l6oi 

3 7. 74 92 

107 .033 

0.92264E-01 

C *4 bub- C i 

2 * 1 343 7 

32. ul 3o 

126.201 

0 . 108 79 

0.450c- Ci 

i • 0660 7 

2 3. Ob 00 

161.618 

0. 1393L 

0.50ut- bl 

i * b b b U 0 

15. 0000 

269.360 

0.23219 

F iL M = 0 . 

i-jOt-bi iN^rtti) 




A, IftLHb 8 

P/P l MI ft) 

PiPSI 

U 1A Vl iFT/ SEb 

MACH NO 

0 

4. 33333 

o 3 . 00 00 

34- . 96 30 

0.301 40E— 0 1 

0 . 5out- b 2 

4.12311 

o 1 . 84 o 6 

36.7476 

0.316 77E-0 i 

b. lUbt- Cl 

3.90137 

38. 3235 

38.8344 

0.334 76E-01 

b. 1 5ub- C 1 

3 • 8668 7 

33.0000 

41 .3223 

0 .356206-01 

b. 2 uut- C l 

3.41303 

3 1. 2348 

44.3391 

0 • 38238E— 01 

0.230b- C i 

3. 1440b 

47.1699 

48 . Id 17 

0. 41533 E- 01 

0- 30bc- C 1 

2 .64duu 

42. 72 00 

53.2005 

0 .458596-01 

L.350t- Cl 

2 . 3 io o 1 

37. 7492 

60. 2060 

0 . 5 1898 E- 0 1 

U.4Ubh- C i 

2.1343 i 

32. 01 3o 

70.9680 

0 .61 1 93E-0 1 

0 .430b- C i 

1 • 60ou 7 

23.0000 

90.9091 

0.78365E-01 

0 . 3bUt- b i 

i.ObOOO 

1 3 * OUOO 

151 .3 15 

0.13061 

F IL M = U . 

iOOb-03 Iftbrltb 




A, iftb He b 

P /P ( Ml ft J 

P,PSi 

U ( A V ) , FT/SEC 

MACH NO 

0 

4. 33333 

63. 0000 

15.5*00 

0.133 96E-0 1 

0. 3U0b- C2 

4.12311 

O i. 8468 

1 6 • 332 4 

0.140 79 E-01 

0. lOub- Cl 

3.90137 

38.5235 

17.2597 

U.14878E-01 

0. l3uc- Cl 

3. 0660 7 

55. OUOO 

id • 3o3 5 

0. L5831E-01 

0 . 2GUE- b 1 

3.41303 

31.2348 

19.7131 

0 • 169 95 E— 01 

b. 2 5bfc- C 1 

3 • i44t»6 

47. 1699 

21.4141 

0 . L8459E-01 

G.3O0b-Cl 

2 . 84800 

42. 7200 

23.6447 

0 .2G3S2E-Q1 

b,3:>Uc- Ci 

2.3 lOol 

37. 7492 

-2 8 . 7 3 8 2 

0 • 23066E— 01 

G.4U0L-C1 

2.13437 

32.0156 

31.5502 

0.27197E-01 

b.430b- b 1 

1 • 0660 7 

25. GOOD 

40.4040 

0 .348 29E-01 

0.300b— C 1 

1 . 000 OU 

15. 0000 

67.3400 

0 .58048E-01 


Ul 


cn 

to 


PLOT Oh POWER VS H(MEAN) 


U.OOUlOO 
0 . 00100 * 

L 

1 

i 

1 

1 

1 

1 

1 

1 

0. 000801 
1 
i 
1 
1 
1 
1 
1 
1 
1 

0. COOoU 1 
1 
l 
l 
1 
1 
1 
1 
1 
i 

0. 000401 
1 
1 
1 
1 
1 
1 
1 
1 
1 

0. 000201 

0. 000100 


0. 0u012 5 0. 000150 0.000175 0.000200 0.000225 0.000250 0.000275 0.000300 0.000325 0.000350 


0.000125 O.OOGlsO 0.JU0175 0.000200 0.000225 0.000250 0.000275 0.000300 0.000325 0.000350 

POWER Urt VERTICAL SCALE - H ( ME A i\ ) UN HORIZONTAL SCALE 

PCWtK IN HORSE POWER - Tu CONVERT TO WATTS, MULTIPLY BY 745.7 
tiMLAM IN INCHES - TO CONVERT TO MEIERS, MULTIPLY BY 2.54E-2 
EUR SHEAR HEAT IN bTU/MIN, MULTIPLY POWER BY 42.42 





CJ1 

CO 


PLUT UP P/PlrtlNJ X Pur HlMcANi = 0.350J0t-J3 

U. G.0050 u.OiOJ 3 .0 150 0.0200 

5. 001 

1 

1 

i 

1 

1 

i 

* 

1 

i * 

4. OUi 

i * 

i 

1 * 

1 

1 

i * 

1 

1 

1 

i. Gui 
1 
1 
1 
1 
1 
i 
1 
i 
1 

2.001 

l 

1 

1 

1 

l 

i 

1 

l 

1 

1.U01 

1 

1 

l 

1 

1 

1 

1 

1 

1 

0. 0.005U 0.0100 0 .0150 0.0 200 

P/PtMiN) UN VERTICAL SCALE — X UN HURIZuNTAL SCALE 
P /P ( M IN ) IS DIMENSIONLESS 

X IN INCHES - TO CONVERT TO METERS, MULTIPLY BY 2.54E-2 


0.0250 


* 


0.0250 


0.0300 0.0350 0.0400 0.0450 


* 


* 


* 


* 


0.0300 0.0350 0.0400 0.3450 


0.3500 


* 


0.0500 



LOMPKESg IBL c SEALING DAM wi TH SMALL f I L F ANGLE 


SAMPLE PrU8i_eM 


TILT ANGLE ^ 0. KADI AN 5 


P2» N/M 2 

PI » N/M2 

T , DEG 

K 

V IgCOS IT Y » N-S /M2 

MOLECULAR WEIGHT 

CP/CV 


u.lu342E 

E6 U • **4 til 6t 06 

3 L 1 

• 

0. 1896 1 £—04 

28.9660 


L. 40000 


K2»Mh TER 

S K i , ME Tt KG 

L, METERS 

KHU, 

K0/M3 

RHU(ROT) ,KG/M3 

NO 

OF GRID POINTS 

0 .8420 lt~ 

G 1 0. 3293 lt-Ol 

0 


0 


0 


n 


N , KP S 

V,M/S 

CP, J/ KG 

-UEG K 


SK IP A 

SKIP R SKIP T 




23.3UGU 

0 

iOO^t. 

7d 


F 

F F 




BEGIN duipu r 

DA Ta 









GAG CUwSTAN f f J/KG-UCD K 

KHCl 11 , KG/ M3 

A(SUUN6 SPEED), M/S 




26 7 

.066 

3 • 03606 


353.591 





L ,ME I ERG 

AREA, M2 


SPEED 

, KPS 

V, M/S 




U .5260o 

0. o6683E- 

•03 

23.3000 

12. 

2 339 




MEAN F ii_M 

M(OUT) 0 

MACrt 

RE IP) 

RE ( R ) 

KNUDSEN 

F STIFF 

xc 

xc 

F 

MET ERs 

KG/ GEC GC MG 

(MAX) 



NUMBER 

NEWTONS KG/M 

METERS 

BAR 

BAR 

C.2 54E-04 






ANALYSIS NOT VALID 




C. 229E-U4 






ANALYSIS NOT VALID 




C • 2G3E— 04 






ANALYSIS NOT VALID 




0. 178E-04 






ANALYSI S NUT VALI D 




0. 152E-G4 






ANALYSIS NOT VALID 




L. 127E-04 






ANALYSIS NUT VALID 




C. i l^E-u** 






ANALYSIS NOT VALID 




0. i02t-04 






ANALYSIS NOT VALID 




0 .889E-05 

u.L3oE-02 u.LilE-02 

0.711 

273.048 

20.889 

0.008 

138.886 -0 

0.431E-03 

0.339 

0.604 

C. 762E-05 

U . d Got- 03 0.o99E-03 

0.322 

171.94 9 

17.905 

0.009 

138.886 0 

0.431E-03 

0.339 

0.604 

U.635E-05 

U.4S5E-C3 G.4G3E-03 

o.3g 3 

99.307 

14.921 

0.011 

138.886 C 

0* 43IE-03 

0.339 

0.6D4 

0.5L6E-O5 

0. 254E-C3 G.207E-03 

0.2 J>2 

50.948 

11.937 

0.013 

138.886 0 

0.431E-03 

0.339 

0.534 

C.381E-GS 

U.iO/fc-G3 0.874E-04 

0.131 

21.494 

8.953 

0.018 

138.886 0 

0. 431E-03 

0.339 

D. 534 

C • 254E— 0 8 

U.317E-04 0.2d9E-04 

0.038 

6.368 

5.968 

0.027 

138.886 -0 

0.431E-03 

0.339 

0.634 

MEAN F ILM , ME JfcK S ' POWER, WATTS 

TOTAL HEAT, 

WATTS 

DfcL( T ) , DEG 

K TORQUE, N-M 




0.254E- 04 






ANALYSIS NOT VALID 




0 . 229E-04 






ANALYSIS NOT VALID 




G •2G3E-U4 






ANALYSIS NOT VALID 




0.178b- 04 






ANALYSIS NOT VALID 




0. 152e-o4 






ANALYSIS NOT VALID 




0 • 127E- 04 






ANALYSIS NOT VALID 




0.U4E-U4 






ANALYSIS NOT VALID 




C • 1C2E- 04 






ANALYSIS NOT VALID 




C.689E-05 

0.21286 


0.21309 


0. 15591 

C.91356E-02 




G.762E-05 

0.24834 


0.24860 


0.26883 

0. 10658E-0I 




0.635E- 05 

0.29800 


0.29832 


0.59893 

0.12790E-01 




G • 5L8E- 05 

0.37261 


0.37290 


1.46222 

0.15987E-QI 




0.381E-G5 

0. 49o67 


0.49721 


4.62135 

0.21 316E-01 




G.254E-G3 

0.74601 


0.74561 


23.3956 

0.3197 5E — 01 






M t AN FILM 


Me rtRb 


MEAN 


M EAN 


0 • tics 3u GE- C5 


X f Mt TtK i 

P/P IM 1NJ 

P.N/M2 


LftAVJ » M/6 tC 

MACH NU 

U 

4. iidJi 

0.4^6 lot 

03 

5 d ,0233 

0.16410 

U.127l-C3 

h. i2di i 

0. 42u42 E 

06 

oO. 98 Id 

0.17246 

0 • 2 5*+t- C5 

3.90157 

0. h 03 3 i E 

Uo 

o4.4444 

0.18226 

0.3ti lc- C3 

8 • 0666 7 

G.i792LE 

Oo 

od, 573U 

0.19393 

O.budt:- Li 

3-4 1 5o5 

U.35325E 

06 

73.6124 

0.20818 

U • o J3 l- L n 

i • N4oo 

G. 320236 

06 

74*9559 

U.220 13 

0 • 7 6<ic- Ci 

2 . 84tiO0 

0.2 44 54 E 

Oo 

dd .28 43 

0.24968 

0.oo9l- Li 

2- 5iooi 

0.2o027t 

06 

99,90^8 

0 .28256 

0.1 02t“ U 

2.13437 

0. 2^o 746 

UO 

117.802 

0.33316 

0*1 l4t- C2 

i • 6060 7 

0. 172 0 7t 

Oo 

1 30 .d6 1 

0.42665 

0. 12 /£- C2 

i . UOOUL 

0. 1 03 42 t 

Oo 

231.434 

0.71109 

F iL M = 0 • 

/o20Ut-C5 

Me I£ft 3 




X, Ml f cK 5 

P/P MINI 

P,N/M2 


UUV J »M/6EC 

MACH NU 

u 

4. iiiii 

0- 448 L 6E 

U6 

42 .o294 

0.12056 

Ci. 12 7t- C3 

+.12311 

0.42642 t 

Gt> 

44 . tiO 29 

0.12671 

o.2:><*l- Ci 

j . 90i 3 / 

0.403516 

Oo 

47.3469 

0.13390 

0. Jdlt-Ci 

i • ouoo 7 

0.3 7421 £ 

0 o 

50 . 38u 1 

0.14248 

0*5 Ude- L 3 

i . 4 L5o3 

0.55323E 

00 

34 .0826 

0.15295 

tl.o t- Ci 

i. I446u 

o. 32523 E 

00 

38.7431 

0.16613 

0* 7o2t- Li 

2 . 64dUU 

0. 244-34E 

Oo 

o4 . 8o2 1 

0 .18344 

C.bd^L- Li 

2 • o iool 

0.2ou27t 

Oo 

73.4031 

0.20759 

U. 1 U2t- C2 

2 . i i4i 7 

6. 220 74 1 

Oo 

do . 548o 

0.24477 

0.1 J.4E- C2 

i .06667 

U. 1 72 3 7t 

06 

1 10 .d 30 

0.31346 

0 • i 27 l- Lc 

i . UGQGU 

0. 1 03 42 t 

06 

184.727 

0.52243 

F ILH = U . 

fci 3U0h— L5 

M l r E K S 




X,Mt Jbrt i 

H/P MINI 

P* N/M2 


OtAV* ,M/S£C 

MACH NO 

0 

4 • 3 i iii 

0. 4*+8 16E 

Oo 

29.0037 

0 .83723E-) 

U. 12/L- Ci 

i2ili 

0.42o42t 

00 

31. 1131 

0.87992E- 

0.2:>4t- Ci 

3. 9C 157 

0. 40331 E 

06 

32 . 87 98 

0.92988E- 

0.38 it- Ci 

i - o<doo7 

0.37921E 

06 

34.9862 

0.98945 E- 

0.5 Udfc- Li 

i • 4 i 5o5 

0. i 3i 2 3 E 

Oo 

37.5374 

0.10622 

0 .o i Ot- C i 

i . 14460 

0. 32323 E 

Oo 

40.7938 

0.11537 

U. 7o2t- Ci 

2.84800 

0. 2 9454E 

Uo 

43.0431 

0.12739 

u. a cJ4l- Ci 

2 . 6 lbol 

U. 2o02 7E 

Oo 

50.9744 

0.144 16 

0 • 10 2t- C2 

2.1343 7 

U.22074E 

06 

O0.1032 

0.16998 

U.l l4t- C2 

1 . 0660 7 

0. 1 72 3 7t 

06 

76.9697 

0.21768 

0.127c-C2 

i • GGuOU 

0. i 03 42 E 

06 

128.283 

0.36280 


cn 

ai 



cn 

05 


M LAN 


MEAN 


MEAN 


F 1l M - 0 

. 50 bUCt- 05 

ME TER 5 




X» M t TEK i» 

P/P (Mi NJ 

P, N/M2 


U(AV) f M/SEC 

MACH NO 

u 

*♦.33333 

0. 448 L b t 

06 

18 .9464 

0 . 5 35 83 E- 01 

a.i^rb-Cj 

4 • i 23 i i 

0. 42642 L 

Ob 

19.9124 

0 • 563 15E-G 1 

0* 2 5**L- C 3 

j.yuib? 

0. 40351 1 

Ob 

21.0431 

0.595 12E-0L 

0. Jo II- 12 

3 . 0 o o o 7 

0. 3792 1 E 

06 

22.3912 

U.63325E-01 

U.5U8C- L 8 

3.4 1 8o5 

0.35325E 

Ob 

24.0367 

0.67979E-01 

U .6 35t- C 3 

3 • l44oo 

0.32523 t 

Ob 

26. 1081 

0. 73837E-01 

0.70 2c— l3 

2 • b4ou0 

0 . 2 94 64 E 

Ob 

28 .827b 

0. 81528 E- 01 

u.ouyt- C j 

2 . blo6i 

0.2602 7 t 

06 

32.623b 

O.92264E-0I 

o. luce- C2 

2.13437 

0. 2 2 0 74 1 

06 

38.466 l 

0.10879 

U.ll4c-C2 

i • b66o 7 

0 . 1 723 7t 

Oo 

49.2b06 

0.13931 

U.127E-C2 

i . OUUOU 

0. 10342 E 

Ob 

82.10 10 

0.23219 

FIlM = U 

. JolOOfc— 05 

Mt IERS 




XfMcFCKS 

P/P (MINI 

P, N/r12 


U(AV), M/SEC 

MACH NO 

0 

4.33333 

0. 448 1 6 1 

06 

10.6573 

0. 30140 E-01 

U. Ic7t- l3 

4.12311 

0. 42b42 t 

Ob 

11.2007 

0 .3 16 77 E-0 1 

U.264C- C3 

:> • 901 ol 

o. 403 5 1 E 

Ob 

11.8367 

0.3 34 76 E-01 

0.3b it- C 3 

3 • ot»bo7 

0.37921 E 

06 

12.5950 

0 .35620 E-01 

0 • 508c- C3 

3 . 4 1 5 o5 

0. 3 53 23 £ 

Ob 

13,5206 

0 • 3B238E-01 

0.6 6 JC- C3 

3 . 14406 

0* 32 323 E 

Ob 

14.6858 

0.41533 E-01 

O.762E-03 

2.848U0 

0.29454c 

Ob 

16.2155 

0.45859E-01 

U.889E- L3 

2 . 5 lb6i 

0. 2 602 7 E 

06 

18.3508 

0.51898 E-01 

0. 102t- C2 

2.1343 7 

0.22074E 

Ob 

21.6372 

0.61193 E-01 

0.U4E-L2 

t • bbfco 7 

0.17237E 

06 

27.7091 

0.78365 E-01 

0.127E-C2 

1.00000 

0.1 0342 E 

06 

46. 1818 

0.13061 

FILM = 0 

• 25400t-05 

METERS 




X» Me TfcKS 

P/P (MINI 

P,N/M2 


lU AVI , M/SEC 

MACH NO 

0 

4.33333 

G.44816E 

06 

4.73660 

0.13396E-0L 

0. 127E-L3 

4.12311 

0. 42642 E 

Ob 

4.978 10 

0.140 79 E-01 

0.254c- 03 

3 • 9 01 d 7 

0.40351 E 

06 

5.26077 

0 • 148 78 E-01 

0.38 It- L3 

3 • 6o6u 7 

0. 37921 E 

Oo 

5.59779 

0 • 15831 E-OL 

0.50BE-CJ 

3.41505 

0.35325 £ 

06 

6 .00918 

0.16995E— 01 

U.635E-C3 

3. 1 44bo 

0.32523E 

06 

6.52701 

0. 184 59 E— 01 

0.7 6 2C- C3 

2.84800 

0. 2 9454 E 

06 

7.2 Qb89 

0 . 20 382 E- 01 

0.889C-C3 

2 . 31bbl 

0.2b027tw06 

8.15591 

0.23066 E-01 

0. 1G2C- C2 

2.1 3437 

0.22074E 

Ob 

9.61651 

0.27197 E-01 

0.1 14E- C2 

1 • 6066 7 

0.17237E 

06 

12.3151 

0 • 348 29 E— 01 

0 • 1 27E- C2 

1 .00000 

0. 1 03 42 E 

06 

20.5252 

0.58048 E-01 



COMPRESS 1 BL t SEALING DAM WITH SMALL TILT ANoLE 


SAMPLE PROBLEM 
INPUT DATA - 


TILT ANGLE = 0.0010 RADIANS 


P2.PSIA 

P I » PSi A 

T.OEG F 

VIS COS I TY ,L 6-SEC/ FT 2 MOLECULAR 

WEIGHT 


CP/CV 

I5.00GC 

65.0000 

100. 

0 

.39600 E-06 

28.9660 


1.40000 

R2, INCHES R I 

.INCHES 

L, INCHES 

RHO ,L B-SEC2/ FT4 

RHOIROT ) , LB- SEC2/F T4 

NO 

OF GRID POINTS 

3 . 3 i 50 C 

3. 2o500 

20.o717 

0 


0 



11 

N , KPM V 

.ft/sec 

C P . 3 T U/ LB— DEG R 


SKIP A 

SKIP R SKIP 

T 



1398. UC 

40. 13 75 

0.24000 


F 

F F 




ifcCiN OUTPUT DATA 









GAS CONSTANT ,FT-LB/Ld (M)-OEG R 

RHGU) . LB-SEC2/ FT4 

AT S OUND SPEED), 

FT/SEC 




53.3522 


0.9737IE-02 


1160 .08 





L , INCHES 

AREA, IN2 

SPEED, RPM 

V, FT/SEC 




20.0717 

1. 03356 

1396.00 


40. 1375 





MEAN FILM M(DUT) 

u 

MACH RE ( P) 

RE ( R) 

KNUDSEN 

F STIFF 

XC 

XC 

F 

INCHtS Lb/MIN 

SCFM 

(MAX) 


NUMBER 

LB LB/IN 

INCHES 

BAR 

BAR 

C.- I0UE-02 





ANALYSIS NOT VALID 




C.9GUE-U3 





ANALYSIS NOT VALID 




C.80GE-G3 





ANALYSIS NOT VALID 




C. 7 COE- 03 





ANALYSIS NOT VALID 




U6U0E-03 





ANALYSIS NOT VALIO 




0.500E-U3 





ANALYSIS NOT VALID 




L.450E-U3 





ANALYSIS NOT VALID 




C. 4G0E- 03 





ANALYSIS NOT VALIO 




C.33DE-U3 0.17/ 

2. 31 o 

0.734 2d 6 . 097 

19.865 

0.008 

29.507 -0.464E 04 

0. 181E-0I 

0.363 

0.571 

0.300E-OJ u.iii 

1.431 

0.33d 162.424 

16.879 

0.009 

29.220 -0.677E 04 

0. L80E-01 

0.361 

0.565 

C.250E-03 0.O3OE-01 

0.632 

0.391 106.207 

1 3.393 

0.Q11 

28.819 -0.954E C4 

0. 1 79E-0I 

0.359 

0.558 

C • 20U t-U 3 0. 320E-C1 

0.41 9 

0.233 34 .07 I 

L0.907 

0.014 

28.221 -0.147E 05 

0. 177E-01 

0.355 

0.546 

C.150E-03 U . i 2 Jt- 0 1 

0.1 70 

0.144 23.044 

7.924 

0.019 

27.240 -0.262E 05 

0. 174E-0I 

0.349 

0.527 

C. IGuE-u3 0 • J 46E— 02 

0.432t-0l 

0.0o4 6.359 

4.960 

0.029 

25.387 -0.506E 05 

0. 169E-01 

0.338 

0.491 

IEAN FILM, INCHtS 

Pu*tR,H. P. 

SHEAR HEAT , 8T U/ M I N 

DtLTT i, DEG F 

TORQUE, FT-LB 




C. I0UE-U2 





ANALYSIS NOT VALID 




t . 90ut- 03 





ANALYSIS NOT VALID 




L.sOOE-03 





ANALYSIS NOT VALID 




C.70OE-O3 





ANALYSIS NOT VALID 




C.6UUE-03 





ANALYSIS NOT VALIO 




U.500E-U3 





ANALYSIS NOT VALID 




G.45UE-U3 





ANALYSIS NOT VALID 




C.4GUE-UJ 





ANALYSIS NOT VALID 




C. 330t-U3 

0. 2o3‘*5t-G3 

0.12L09E-OL 


0.2849 7 

0.6738 IE-02 
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